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IMPROVED  DRAINAGE  AND  FROST  ACTION  CRITERIA 
FOR  NEW  JERSEY  PAVEMENT  DESIGN 
PHASE  2:  FROST  ACTION 


by 

Richard  L.  Berg  and  Richard  W.  McGaw 


INTRODUCTION 

In  early  1975  engineers  from  the  New  Jersey  Department  of  Trans- 
portation (NJDOT)  visited  CRREL  to  discuss  a program  they  were  embarking 
upon  with  the  assistance  of  the  Federal  Highway  Administration.  The 
NJDOT  was  interested  in  including  open-graded  drainage  layers  in  their 
pavements  to  remove  excess  water.  Before  constructing  actual  field 
installations,  they  chose  to  determine,  analytically,  the  effect  of 
frost  penetration  beneath  pavements  when  using  an  open-graded  drainage 
layer.  In  early  July  1975,  the  NJDOT  issued  a contract  to  CRREL  to 
conduct  this  study,  the  principal  goal  of  which  is  to  determine  analyt- 
ically the  impact  of  an  open-graded  drainage  layer  on  the  depth  of  frost 
penetration  beneath  pavements  in  New  Jersey.  The  actual  effect  of  a 
particular  drainage  layer  will  depend  on  its  thickness  and  location  in 
the  pavement  profile,  the  gradation  of  the  soil  and  the  stabilizing 
agent  applied,  the  range  of  water  contents  likely  to  be  encountered,  and 
the  resulting  thermal  properties. 

The  study  had  the  following  goals: 

a.  To  determine  whether  the  depth  of  frost  penetration  would 
change  significantly  if  an  open-graded  drainage  layer  was  used,  and 

b.  To  determine  whether  exceptional  adverse  effects  should  be 
guarded  against. 

The  modified  Berggren  equation  was  used  to  estimate  the  maximum 
seasonal  frost  penetration  depths  beneath  several  pavement  profiles. 
Previous  experience  has  indicated  that  this  equation  predicts  frost 
depths  which  are  generally  adequate  for  engineering  applications.  To 
provide  a higher  degree  of  reliability  in  the  predictions,  however,  it 
was  necessary  to  measure  surface  transfer  coefficients  for  selected  New 
Jersey  pavements  and  to  conduct  thermal  conductivity  tests  on  repre- 
sentative base  and  subbase  course  materials,  pavement  materials  and 
open-graded  drainage  materials. 

Air  and  surface  temperature  data  from  the  1975-1976  winter  were 
supplied  by  the  NJDOT  and  were  used  to  determine  surface  transfer 
coefficients.  The  NJDOT  also  furnished  aggregate  for  the  thermal  con- 
. ductivity  tests  on  the  base  and  subbase  course  materials  and  specimens 

for  thermal  conductivity  tests  on  the  paving  materials . 
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SURFACE  TRANSFER  COEFFICIENTS 


To  calculate  seasonal  frost  penetration  depths,  it  is  necessary  to 
determine  the  temperature  condition  at  the  ground  surface.  Since  air 
temperatures  are  generally  available  and  surface  temperatures  are  not , a 
correlation  between  these  temperatures  is  required  to  establish  the 
thermal  boundary  condition  at  the  ground  surface.  Unfortunately  no 
simple  correlation  exists  between  air  and  surface  temperatures.  The 
difference  between  air  and.  surface  temperatures  at  any  specific  time  is 
influenced  by  latitude,  cloud  cover,  time  of  year,  time  of  day,  atmospheric 
conditions,  wind  speed,  exposure,  surface  characteristics  and  subsurface 
thermal  properties. 

Air  and  surface  temperatures  are  converted  to  "freezing  indexes"  to 
compute  seasonal  frost  penetration  depths.  To  understand  the  concept  of 
freezing  indexes,  two  definitions  are  necessary: 

(l)  Degree-days . The  degree-days  for  any  one  day  equal  the 
difference  between  the  average  daily  temperature  and  32°F.  Freezing 
degree-days  occur  when  the  average  daily  temperature  is  below  32°F. 

Thawing  degree-days  occur  when  the  average  daily  temperature  is  above 
32  °F. 


(2)  Freezing  index.  The  freezing  index  is  the  number  of  degree- 
days  between  the  highest  and  lowest  points  on  a curve  of  cumulative 
degree-days  versus  time  for  a freezing  season.  It  is  a measure  of  the 
duration  and  magnitude  of  below- freezing  temperatures  for  a given 
freezing  season.  Air  temperatures  for  the  air  freezing  index  are 
measured  approximately  1*.5  ft  above  the  ground  while  temperatures  for 
the  surface  freezing  index  are  measured  immediately  below  the  surface. 

The  ratio  of  surface  freezing  index  to  air  freezing  index  is 
designated  as  the  "n- factor"  for  freezing  conditions.  Determination  of 
the  n-factor  for  a specific  location  requires  concurrent  observations  of 
air  and  surface  temperatures  throughout  several  complete  freezing 
seasons. 

Lunardini  (1977)  has  conducted  an  extensive  study  of  n-factors . 

For  asphaltic  concrete  pavements  under  freezing  conditions,  he  found 
that  the  range  of  values  from  25  different  sites  was  from  1.00  to  0.25 
with  an  average  of  0.57.  For  Portland  cement  concrete  pavements  sub- 
jected to  freezing  conditions,  the  range  of  values  was  0.95  to  0.12  for 
21  sites  and  the  average  value  was  0.58.  Moulton  (1968)  studied  frost- 
related  problems  in  West  Virginia  pavements.  He  reported  values  of  O.96 
to  0.25  for  asphaltic  concrete  pavements  in  the  winter  and  O.87  to  0.12 
for  Portland  cement,  concrete  pavements.  It  is  the  authors'  opinion  that 
a substantial  amount  of  the  variation  in  n-factors  is  due  to  differences 
in  placement  of  the  surface  temperature  probe  and  differences  in  proximity 
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of  the  surface  and  air  temperature  observation  locations.  Recent  studies 
at  CRREL  by  Eaton  (1976)  yielded  n-factors  of  0.57  to  0.1*3  for  asphaltic 
concrete  pavements.  The  Department  of  the  Army  (1966)  recommended 
design  values  of  0.9  for  "both  asphaltic  concrete  and  Portland  cement 
concrete  pavements  kept  free  of  ice  and  snow.  These  values  have  been 
used  for  many  years  for  computing  maximum  seasonal  frost  depths  for 
pavement  design.  The  larger  values  for  n-factors  would  result  in  greater 
calculated  frost  penetration  depths  because  the  surface  freezing  indexes 
would  be  larger. 

Since  the  range  in  values  for  n-factors  is  quite  large,  the  CRREL 
investigators  suggested  that  pavement  surface  temperatures  be  measured 
to  obtain  representative  data  for  New  Jersey.  During  the  1975-1976  and 
1976-1977  freezing  seasons,  the  New  Jersey  DOT  measured  air  and  surface 
freezing  indexes  on  Portland  cement  concrete  and  asphaltic  concrete 
pavements  at  three  different  locations  (Figure  l).  Site  1 was  located 
on  Interstate  295  near  Bordentown,  Site  2 was  on  Interstate  287  near 
Bedminster,  and  Site  3 was  on  Route  15  near  Rockaway.  Figures  2 - U 
contain  air  and  surface  freezing  indexes  from  Sites  1,  2 and  3,  respec- 
tively for  the  1975-1976  winter.  Data  for  the  1976-1977  winter  are  in 
Berg  (in  preparation).  Both  Portland  cement  concrete  and  asphaltic 
concrete  pavement  surface  indexes  are  shown.  Figure  2 illustrates  the 
method,  discussed  above,  for  determining  the  air  freezing  index.  Other 
freezing  indexes  are  determined  using  the  same  procedure. 

Three  days  of  data  are  missing  in  mid-December  at  the  Rockaway  site 
and  approximately  15  days  of  data  are  missing  in  late  November  and  early 
December  at  the  3edminster  site.  The  freezing  season  probably  began 
during  the  time  of  missing  data  at  both  of  these  sites;  however,  n- 
factors  determined  from  the  partial  freezing  seasons  are  sufficiently 
accurate  for  the  present  studies. 

Table  1 contains  freezing  index  values,  n-factors  and  the  length  of 
the  freezing  season  measured  at  each  of  the  three  sites  during  the  1975- 
1976  winter.  Bordentown  was  the  warmest  site,  with  an  air  freezing 

Table  1.  Freezing  indexes  and  n-factors 
1975-1976  winter 

Air  Asphaltic  Concrete  Portland  Cement  Concrete 


Site 

Season 

Index 

Season 

Index 

n-factor 

Season 

Index 

n-factor 

Days 

°F  Days 

Days 

°F  Days 

Days 

°F  Days 

Bordentown 

lo 

95 

20 

59 

0.62 

20 

57 

0.60 

Bedminster 

1*0 

151 

1*0 

127 

0.8U 

Uo 

82 

0.51* 

Rockaway 

58 

1*71 

56 

216 

0.1*6 

56 

188 

0.1*0 

index  of  95  °F-days,  and  Bedminster  and  Rockaway  had  air  freezing  indexes 
of  151  °F-days  and  1+71  °F-days , respectively.  At  all  three  sites  the  n- 
factors  for  the  asphaltic  concrete  surface  were  larger  than  those  of  the 
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Figure  U.  Air  and  pavement  surface  freezing  indexes  near  Rock- 
away,  1975-1976  winter. 

Portland  cement  concrete  surfaces.  The  differences  were  rather  small, 
however,  at  the  Bordentown  and  Rockaway  sites.  The  n-factor  for  the 
asphaltic  concrete  pavement  at  Bedminster  is  significantly  larger  than 
the  remaining  data.  Data  from  Figure  3 show  no  anomalous  behavior  of 
the  asphaltic  concrete  temperature  sensor  at  Bedminster,  and  NJDOT 
engineers  (Cosaboom,  1977)  state  that  the  instrumentation  functioned 
properly  during  the  winter.  The  air  freezing  index  measured  at  Bedminster 
during  the  1976-1977  winter  was  UU6°F-days,  and  n-factors,  (Berg,  in 
prep.)  were  0.90  and  0.87,  respectively,  for  the  asphaltic  and  Portland 
cement  concrete  pavements.  Air  freezing  indexes  measured  at  National 
Oceanic  and  Atmospheric  Administration  (formerly  U.S.  Weather  Service) 
sites  at  Long  Valley  and  Sommerville,  which  are  approximately  9 miles 
north  and  9 miles  southeast  of  Bedminster,  respectively,  were  800  °F- 
days  and  790  °F-days.  If  the  average  of  795  °F-days  is  used  with  the 
measured  surface  temperatures  at  Bedminster,  the  1976-1977  n-factors 
reduce  to  0.50  and  0.U9  for  the  asphaltic  and  Portland  cement  pavements. 
These  values  agree  much  more  closely  with  n-factors  from  the  other  two 
sites  in  New  Jersey.  The  n-factors  used  in  calculating  frost  depths  in 
this  study  are  shown  subsequently.  The  effect  of  different  n-factors 
is  discussed  further  under  the  section  "Calculated  Frost  Depths". 
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THERMAL  CONDUCTIVITY  TESTS 


Materials  Tested 


The  NJDOT  shipped  samples  of  typical  base  course,  pavement,  and 
candidate  drainage-layer  materials  to  CRREL  for  thermal  conductivity 
tests.  The  base  course  soils  were  subjected  to  the  Guarded  Hot  Plate 
(GHP)  test,  while  the  Probe  test  was  used  for  the  pavement  and  stabilized 
drainage-layer  materials.  One  drainage-layer  material  was  tested  in  the 
unstabilized  condition,  using  both  the  GHP  and  the  Probe  tests. 

Appendix  A contains  data  sheets  describing  in  detail  the  physical 
properties  of  all  the  materials  as  they  were  received  at  CRREL.  Several 
of  these  properties  are  also  discussed  in  the  following  sections. 

Base  Course  Soils  (Guarded  Hot  Plate  Test) 


Four  base  course  materials  were  received  from  the  NJDOT  in  bulk, 
and  were  compacted  at  CRREL  to  approximately  95%  of  the  maximum  density 
at  optimum  moisture  content  listed  in  Table  2.  The  maximum  dry  density 
and  optimum  moisture  content  values  were  obtained  from  data  provided  by 
the  NJDOT  (see  Fig.  Al).  Each  material  was  tested  at  two  water  contents 
in  both  the  frozen  and  thawed  states.  Sample  size  for  the  Guarded  Hot 
Plate  apparatus  was  20  in.  x 20  in.  x 3 in.  thick.  A complete  descrip- 
tion of  the  test  equipment  and  procedure  is  provided  in  Appendix  B. 
Physical  properties  of  the  base  course  materials  are  summarized  in  Table 
2. 


Table  2.  Fhysical  properties  of  base  course  materials. 


Actual  % 

Passing 

Sieve 

Ogdensburg 

Hamilton 

Hamburg 

Pennington 

Size 

Gneiss 

Lakes 

Shale 

Traprock 

4" 

_ — _ 

100 

_ __ 

2" 

100 

— 

100 

100 

3/4" 

86.9 

98.4 

93.2 

83.1 

tf  4 

48.0 

74.4 

69.1 

55.1 

K 50 

9.5 

11.9 

7.2 

16.8 

0200 

3.8 

1.4 

2.8 

8.9 

Opt inum 
Moisture  (%) 

8.1 

7.6 

9.5 

11.1 

Max . Dry 
Density  (pcf) 

135.6 

123.6 

127.8 

135.9 

Minimum 

Porosity 

0.210 

0.225 

0.225 

0.244 

The  two  moisture  contents  chosen  for  the  materials  were  (l)  approxi- 
mately 80#  of  the  optimum  moisture  content,  and  (2)  the  "field  capacity" 
of  each  material..  The  former  was  intended  to  simulate  a rather  dry  i 

field  condition,  and  the  latter  a nearly  saturated  field  condition. 

Initial  plans  were  to  conduct  the  higher  moisture  content  tests  at  full 
water  saturation;  however,  it  was  not  possible  to  attain  saturation  , 

because  water  tended  to  flow  out  of  the  mold.  Field  capacities  were 
obtained  by  wetting  the  3-inch  thickness  of  soil  with  as  much  water  as 
it  would  hold  without  draining. 

The  two  temperatures  chosen  for  each  set  of  thermal  conductivity 
tests  were  approximately  UO°F  and  22°F.  These  temperatures  were  chosen 
to  provide  a value  in  the  frozen  state  and  one  in  the  thawed  state , and 
to  faciliate  a comparison  with  the  values  reported  by  Kersten  (19^9). 

Results  of  the  Guarded  Hot  Plate  tests  (Table  3 and  Fig.  5)  indicate 
that  three  of  the  four  unstubilized  New  Jersey  soils  have  higher  conduc- 
tivities in  the  thawed  state  vhan  in  the  frozen  state.  Kersten  (19^9) 
observed  a similar  behavior  in  soils  he  tested.  Discussing  differences 
in  thermal  conductivity  at  25°F  and  at  UO°F  (Fig.  6)  Kersten  stated: 

"It  was  found  that  the  difference  in  K values  at  these 
two  points  was  chiefly  dependent  upon  the  moisture  content. 

On  relatively  dry  soils,  for  example  those  in  the  air-dry 
condition,  very  little  difference  was  found.  As  the  moisture 
content  was  increased,  the  K at  25  degrees  became  less  than 
that  above  freezing.  With  a further  increase  in  moisture 
content,  the  below- freezing  value  became  progressively  greater 
than  that  above  freezing.  The  ratio  of  the  frozen  to  unfrozen 
value  also  depended  somewhat  upon  the  density  of  the  soil,  the 
ratio  in  general  being  greater  for  a high  density  test  than  for 
a low  density  test  at  the  same  moisture  content." 

Probably  the  Ogdensburg  Gneiss,  Hamilton  Lakes  Sand  and  Hamburg 
Shale  were  tested  at  moisture  contents  where  the  ratio  of  the  thermal 
conductivity  in  the  frozen  state  to  that  in  the  thawed  condition  was 
less  than  unity,  but  the  Pennington  Traprock  was  tested  at  moisture 
contents  where  the  ratio  was  greater  than  one.  Tests  at  a series  of 
moisture  contents  would  be  required  to  develop  the  relationships  for 
these  conditions  but  were  unnecessary  for  this  study. 

The  results  of  the  GHP  tests  are  compared  with  Kersten 's  mean 
values  for  sandy  soils  in  a later  section.  They  are  also  used,  along 
with  the  results  of  the  Probe  tests,  in  the  computation  of  frost  pene- 
tration which  will  be  discussed  subsequently. 
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Pavement  Materials  (Probe  Tests) 


Two  specimens  each  of  Portland  cement  concrete  and  of  three  types 
of  bituminous- stabilized  pavement  materials  were  received  from  the  NJDOT 
and  tested  by  the  Probe  method.  A detailed  description  of  the  probe 
theory,  sample  preparation  and  experimental  techniques  is  presented  in 
Appendix  C. 

The  bituminous  mixtures  included  an  asphaltic  surface  course 
(Mix  #5),  a binder  course  (Mix  #2),  and  a base  course  (Mix  #l).  All  the 
asphaltic  concrete  mixes  contained  about  3%  bitumin  (AC-20)  by  weight. 

All  of  the  materials  were  tested  in  the  as-received  condition,  (air- 
dry)  both  above  and  below  freezing.  Normally  one  sample  of  each  material 
was  used  for  the  Probe  tests.  Properties  of  the  pavement  materials  are 
given  in  Table  U.  Properties  of  the  stabilized  open-graded  drainage 
layer  materials  are  also  shown  in  the  table  but  are  discussed  in  the 
next  section. 

The  Portland  cement  concrete  samples  had  been  prepared  from  a 
paving  mixture  on-site.  The  bituminous  concrete  samples  were  prepared 
in  the  NJDOT  soils  laboratory.  Samples  were  cylinders  6 in.  in  diameter 
and  approximately  8 in.  long.  A steel  rod  0.0T0  in.  in  diameter  was 
molded  along  the  axis  of  each  cylinder;  it  was  removed  prior  to  testing 
and  the  thermal  conductivity  probe  inserted  in  its  place. 

Prior  to  testing,  the  specimens  were  brought  to  the  test  temperature 
with  the  probe  in  place.  A pre-selected  heating  current  was  applied  for 
a period  of  10  or  15  minutes,  during  which  the  temperature  of  the  probe 
was  recorded.  The  current  was  then  turned  off  and  the  specimen  allowed 
to  cool.  Results  obtained  during  both  heating  and  cooling  were  averaged 
to  obtain  the  value  of  thermal  conductivity  used  in  the  pavement  analysis; 
normally  heating  and  cooling  values  were  within  3%  of  each  other.  Table 
5 presents  the  thermal  conductivity  test  results  for  the  pavement  materials . 
These  data  are  plotted  in  Figure  7- 

Results  from  the  tests  indicated  that  the  bituminous  base  course 
(Mix  # l)  and  the  binder  course  (Mix  #2)  had  essentially  the  same  thermal 
properties;  therefore  they  were  thereafter  treated  as  one  material.  The 
results  also  indicate  that  the  thermal  conductivity  of  New  Jersey  pave- 
ment materials  increases  with  lower  temperatures,  including  temperatures 
below  the  freezing  point  of  water.  This  trend  is  similar  to  results  from 
Pennington  traprock  and  may  be  the  result  of  similar  behavior  except  that 
the  water  is  replaced  by  cement  mortar  or  asphalt  cement. 

Bituminous-stabilized  Ooen-graded  Drainage  Layer  Materials  (Probe  Test) 

Three  gradations  of  bituminous-stabilized  open-graded  drainage-layer 
materials  (crushed  Mt.  Hope  granite)  were  prepared  by  the  NJDOT.  Mixes 
ttl  and  H2  contained  U*  bitumin  (AC-20)  by  weight,  and  Mix  #3  contained 
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Table  Physical  properties  of  New  Jersey 

pavement  materials  for  probe  tests . 


U.S.  Sieve 
Size 

Percentage  Finer  by  Weight 

P.C.  concrete 
mixture 

Coarse  Fine  Comb'd, 

Bit.  concrete 
mixtures 

#1  #2 
Base  Binder 

#5 

Surface 

Open-graded 
bit.  mixtures 

#1  #2  (not 

tested) 

#3 

mm 

% 

Of 

/0 

% 

% 

/0 

(0 

% 

% 

% 

% 

2"  50.8 

100 

1-1/2  38.1 

100 

100 

100 

1"  25. 4 

100 

100 

99.8 

99.2 

- 

- 

3 A"  19.0 

- 

95.2 

- 

80 

95 

1/2"  12.7 

32 

56 

79.0 

70.1 

100 

- 

- 

100 

3/8" 

9.5 

20 

100 

49 

70.7 

- 

98.6 

50 

65 

70 

# 4 

4.76 

4 

99 

38 

39.2 

33.0 

64.5 

25 

35 

15 

#8 

2.38 

1 

92 

34 

28.8 

16.9 

50.5 

- 

- 

5 

#l6 

1.19 

- 

78 

28 

- 

- 

- 

- 

- 

- 

mo 

0.595 

- 

49 

18 

- 

- 

- 

- 

- 

- 

mo 

0.297 

- 

16 

6 

11.3 

3.1 

17.2 

3 

6 

0 

#100 

0.149 

- 

4 

1.4 

- 

- 

- 

- 

- 

- 

#200 

0.074 

- 

1 

0.4 

5.2 

1.0 

5.7 

1 

2 

Specific 

Agg. 

2.90 

2.65 

! 

2.81 

2.83 

2.84 

2.76 

1 

2.68 

2.68 

2.6; 

Gravity 

Bind* 

»r 

P. 

I 

C.  = 3.15 

AC-20  = 1. 

,01 

AC-20  = 

1.01 

'omb'  i 

a Solids* 

2.87 

2.60 

2.60 

2.54 

2.55 

2.55 

2. 5' 

Design 

P. 

C.  = 18  JS 

Agg.  = 95# 

Agg.  : 

= 96# 

91% 

Mixture 

Fine  Agg.  = 

29$ 

AC-20  = 5% 

AC-20 

= h% 

3% 

(by  vt ) 

Coarse 

-Ags- 

= 53#  ! 

^Calculated  for  design  mixture. 
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Table  5.  Thermal  conductivity  of  New  .Jersey 
pavement  materials,  probe  method. 


Material 

Cement/ 

Asphalt 

Content 

Aggregate 

Voids 

(Total) 

Air 

Voids 

Aggre- 

gate 

Density 

Overall 

Density 

Heat 

Input 

Thermal 

Conductivity 
(Note  1) 

jf 

~r~ 

ocf 

pcf 

°F 

Btu 

3tu 

mm 

ft -hr 

hr-ft-4? 

r.C.  Concrete 

18.0 

32.0 

119-0 

147-3 

51-3 

0.28 

0.59 

(Surface  Course) 

50.7 

0.29 

0.62 

Air- dry 

i 

50.7 

0.29 

0.60 

50.7 

0.18 

0.57 

(Water /Cement ) 

Patio  6.0'- 

23.4 

0.58 

0.91 

23.4 

0.58 

0.87 

SI 

23.4 

0.58 

0.85 

Bit.  Concrete 

5.0 

18.2 

6.0 

144.8 

152.3 

50.5 

m 

0.55 

#1  (Base) 

50.2 

tm 

0.55 

02  (3inder) 

Air-dry 

23-9 

0.29 

0.66 

23.9 

0.09 

0.65 

3it.  Concrete 

5.0 

17.8 

6.0 

141.5 

149.0 

50.9 

0.55 

0 5 (Surface) 

50.5 

0.55 

Air-dry 

50.5 

0.55 

23-5 

0.22 

0.74 

23.4 

0.31 

0.72 

Open-Graded 

u.o 

28.2 

20.4 

120.0 

124.9 

37.8 

0.30 

0.69 

Bit.  Mix  01 

(Sample 

37.0 

0.30 

0.68 

Drainage ) 

A) 

Air-dry 

24.4 

0.29 

0.73 

24.4 

0.29 

0.71 

4.0 

26.3 

l8.4 

123.1 

128.0 

23.9 

0.29 

0.73 

(Sample 

23.9 

0.43 

0.75 

B) 

Open-graded 

3.0 

30.7 

25.0 

115.9 

119.5 

38.1 

0.15 

0.47 

Bit.  Mix  03 

(Drainage) 

23.4 

0.10 

0.52 

Air-dry 

23.4 

0.15 

0.53 

23.4 

0.15 

0.45 

23.4 

0.22 

0.48 

1 

23.4 

0.22 

0.47 

NOTE  1:  Average  of  heating  and  cooling 

values . 
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Figure  7.  Thermal  conductivity  of  New  Jersey  pavement 
materials , Probe  method . 

3*  bitumen.  Mix  #2  was  eventually  eliminated  from  the  testing  program 
because  the  NJDOT  felt  that  its  density  was  too  high  and  the  resulting 
water  permeability  was  too  low  for  a suitable  free  draining  material. 

Two  samples  of  open-graded  Mix  #1  (OG-l)  and  one  sample  of  open- 
graded  Mix  # 3 (OG-3)  were  tested,  above  and  below  freezing,  by  the  Probe 
method.  Tests  were  performed  in  the  same  manner  as  described  for  the 
pavement  mixtures.  Samnles  were  6-in. -diameter  cylinders  8 in.  long. 

All  tests  were  in  the  as-received  air-dry  condition. 

Physical  properties  of  the  bituminous  stabilized  open-graded  mix- 
tures are  given  in  Table  4;  test  results  are  in  Table  5.  These  results 
are  shown  graphically  in  Figure  7 to  provide  a comparison  with  the  pave- 
ment mixtures. 

It  should  be  noted  that  in  the  above-freezing  tests , the  pavement 
materials  were  tested  near  50°F  to  allow  for  higher  temperatures  caused 
by  solar  radiation  absorbed  at  the  surface.  Because  the  open-graded 
mixtures  would  always  be  at  a lower  depth  in  the  profile,  they  were 
tested  at  a temperature  near  UO°F. 

The  test  results  indicate  that  the  thermal  conductivity  of  OG-3  is 
about  one-third  lower  than  that  of  OG-l,  both  above  and  below  freezing. 
Since  the  aggregate  material  is  the  same  in  both  mixtures,  the  difference 
is  evidently  the  combined  effect  of  a lower  dry  density  and  lower  asphalt 
content  (115.9  pef  and  3 % asphalt)  for  OG-3  as  compared  with  OG-l  (121.5 
pcf  and  b%  asphalt).  These  parameters,  in  turn,  are  controlled  by  the 
differences  in  gradation. 
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On  the  other  hand,  the  thermal  conductivity  of  OG-1  is  similar  to 
the  pavement  mixtures,  which  have  much  higher  dry  densities  (lU2  to  1U5 
pcf)  and  asphalt  contents  of  5$ • In  this  case,  the  higher  thermal  con- 
ductivity of  the  quartzitic  Mt.  Hope  aggregate  evidently  is  sufficient 
to  overcome  the  deficits  in  thermal  conductivity  associated  with  lower 
density,  i.e.  more  voids,  and  a somewhat  lower  asphalt  content. 

Unstahilized  Open-graded  Mix  #1  (Probe  and  Guarded  Hot  Plate  Tests ) 

An  aggregate  having  the  same  gradation  as  that  used  in  the  bitu- 
minous-stabilized open-graded  Mix  #1  (OG-l)  was  tested  in  a moist  con- 
dition. Its  use  in  this  unstabilized  condition  would  be  as  a drainage 
layer  placed  immediately  beneath  a rigid  pavement . 

To  provide  a comparison  between  test  results  from  the  Probe  method 
and  the  Guarded  Hot  Plate  method,  the  unstabilized  OG-1  was  tested  using 
both  techniques.  A slight  difference  in  gradation  occurred  because  of 
the  difference  in  the  size  of  specimen  used  for  each  of  the  tests:  for 
the  Probe  test  specimen,  material  between  the  2-in.  and  the  3A-in. 
sizes  was  replaced  with  3A-in.  material.  The  amount  passing  each  of 
the  smaller  sieve  sizes  normally  would  have  remained  unchanged  by  this 
procedure;  in  actuality,  the  probe  test  sample  was  generally  finer  than 
the  GHP  sample  (Appendix  A).  Apparently  some  breaking  of  corners 
occurred  during  previous  handling. 

This  series  of  tests  was  performed  both  above  and  below  freezing, 
by  each  test  method.  Dry  density  was  121.2  pcf  (GHP)  and  122.5  pcf 
(Probe);  water  content  was  the  "field  capacity"  as  defined  above,  which 
was  approximately  U.5$  on  a dry  weight  basis  f'-r  both  specimens. 

Results  obtained  for  the  unstabilized  open-graded  material  are 
summarized  in  Table  6.  For  comparison,  results  found  by  Kersten  for 
crushed  granite  having  a density  of  120  pcf  and  a moisture  content  of 
and  a sandy  soil  at  a moisture  content  of  i*.5$,  are  also  listed. 

These  data  are  plotted  in  Figure  8. 

The  GPH  values  at  22°F  and  1*0°F  are  between  70$  and  75$  of  the  Probe 
values.  The  reason  is  due  partially  to  the  somewhat  finer  gradation  of 
the  Probe  sample,  and  also  to  the  skin  effect  in  the  GHP  sample,  which 
tends  to  reduce  the  effective  heat  flow  across  the  sample.  The  slight 
difference  in  density  should  have  little  effect. 

In  a comparison  between  the  stabilized  and  unstabilized  drainage 
layer  material  (both  tested  by  the  Probe  method  and  at  the  same  density), 
the  frozen  OG-1  which  was  stabilized  had  a 50$  greater  thermal  conduc- 
tivity than  the  frozen  OG-1  which  was  unstabilized.  In  contrast,  the 
unfrozen  mixtures  had  such  similar  thermal  conductivities  that  two  other 
test  values  were  obtained  (at  approximately  50°F  and  65°F).  The  similarity 
persisted  throughout  the  unfrozen  zone;  thus,  water  and  asphalt  seem  to 
have  nearly  the  same  effect  on  the  thermal  conductivity  of  this  open-graded 
aggregate. 
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Table  6.  Thermal  Conductivity  of  Unstabilized 
Open-Graded  Mix  til. 
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NOTE  1:  Average  of  heating  and  cooling  values. 


W/m°K 


Kersfen 


Crushed  Granite 
«*  = 4.0  % 


Bit.  0G  I 
asphalt*  4.0% 


Probe 


Unstab.  OG  I 
w • 4.5  % 


Unstab.  OG  I 
w * 4.5  % 


Temperature 

Thermal  conductivity  of  unstahilized  open-graded  Mixture 
#1  (w  = water  content)  by  two  methods,  compared  with  the 
corresoonding  bituminous  mixture  and  two  Kersten  materials 


Comparison  with  Kersten's  Values 


Kersten  (19^9)  has  given  equations  and  graphs  for  estimating  the 
thermal  conductivity  of  two  main  groups  of  soils:  sandy  soils,  and  silt 
and  clay  soils.  He  states  that  the  accuracy  is  approximately  ±25%  for 
an  unknown  soil  of  similar  gradation,  density,  and  water  content. 


The  New  Jersey  base  course  aggregates  tend  to  be  sandy  in  character 
and  are  therefore  compared  here  with  Kersten's  predicted  values  for 
sandy  soils.  Gradation  specifications  are  given  in  Appendix  A,  and 
gradation  curves  are  plotted  in  Figure  9. 


US  Std 

Sieve  Opening  U S Std  Sieve  No 


e Ogdensburg  Gneiss  (IA) 

♦ Homllton  Lakes  (1C) 

♦ Hamburg  Shot*  (1C) 

♦ Pennington  Traprock  (5A) 

Figure  9*  Gradations  of  New  Jersey  base  course  materials  used 
in  the  thermal  conductivity  tests. 

The  ratios  of  measured  thermal  conductivities  (K)  to  Kersten’s 
values  (K  ) are  plotted  in  Figure  10  for  the  guarded  hot  plate  tests  on 
the  unfrozen  soils  (Uo°F).  Similar  ratios  for  the  frozen  soils  (22°F) 
are  plotted  in  Figure  11.  For  convenience,  the  ratios  are  plotted 
against  the  percentage  passing  the  No.  200  sieve,  resulting  in  a rela- 
tionship which  should  be  typical  of  New  Jersey  soils.  Table  7 gives 
supporting  data  for  the  two  graphs. 

The  Hamilton  Lakes  sand  (a  well-graded  sand)  is  the  only  one  of  the 
New  Jersey  soils  which  had  thermal  conductivity  values,  both  frozen  and 
unfrozen,  similar  to  Kersten' s data.  The  other  three  soils  have  con- 
ductivities ranging  from  5055  to  75$  of  the  corresponding  Kersten  values. 
These  results  are  nearly  as  one  would  expect,  inasmuch  as  Kersten  states 
that  his  relations  are  for  soil  aggregates  high  in  quartz  content.  The 
Hamilton  Lakes  Sand  is  quartzitic,  whereas  the  Hamburg  Shale,  Ogdensburg 
Gneiss,  and  Pennington  Traprock  have  low  quartz  contents.  As  Kersten 
showed,  a quartz  aggregate  has  the  highest  thermal  conductivity;  a trap- 
rock aggregate  has  the  lowest. 

The  four  New  Jersey  soils  considered  in  these  tests  were  chosen  to 
be  representative  of  the  majority  of  New  Jersey  soils  used  as  base  and 
subbase  course  materials.  Consequently,  it  is  suggested  that  the  following 
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(%)  Passing  No.  200  Sieve  (NJ  Materials) 


Figure  10.  Comparison  of  thermal  conductivity  (K)  of  New  Jersey  soils 

with  Kersten's  values  (K  ) for  sandv  soils,  unfrozen. 

s 


2 4 6 8 10  20 

(%)  Passing  No  200  Sieve  (NJ  Materials) 


Figure  11.  Comparison  of  thermal  conductivity  (K)  of  New  Jersey  soils 
with  Kersten's  values  ( K ) for  sandy  soils,  frozen. 
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procedure  be  used  to  estimate  the  thermal  conductivity  of  other  New 
Jersey  base  and  subbase  course  soils : 


a.  For  the  new  soil,  measure  or  estimate  the  dry  density  and 
water  content. 

b.  Using  the  equations  below,  or  Figures  12  and  13,  find  Kersten's 
value  of  thermal  conductivity  (K  ) for  a soil  having  the  same 
density  and  water  content. 

For  unfrozen  sandy  soils : 

K = (0.7  log  w + O.U)  10°'0lY 
s 

For  frozen  sandy  roils : 

K = 0.076  (10)°'°13y  + 0.032  (10)°-011*6y  (w) 
s 

2 

where  K = thermal  conductivity,  Btu-in/ft  hr  °F 

s q 

Y = dry  density  of  soil,  lb/ft 
w * moisture  content,  $ dry  wt. 

c.  Enter  the  appropriate  curve  in  Figure  10  or  Figure  11  (unfrozen 
or  frozen)  with  the  percentage  passing  the  No.  200  sieve  for 
the  New  Jersey  soil;  find  the  corresponding  ratio  (K/Ks). 

d.  Multiply  the  Kersten  value  (K  ) by  the  ratio  determined  in 
step  (c);  the  resulting  valueSsnould  be  within  approximately 
+10$  of  the  thermal  conductivity  of  the  New  Jersey  soil. 

Kersten  measured  the  thermal  conductivity  of  only  one  bituminous 
paving  mixture.  It  had  a dry  density  of  138  pcf  and  an  asphalt  content 
of  about  6$.  The  measured  values  for  the  New  Jersey  paving  materials 
are  approximately  70$  of  Kersten's  value.  This  value  is  reasonable, 
considering  that  the  New  Jersey  materials  have  asphalt  contents  of  about 
5$,  and  the  aggregates  have  lower  thermal  conductivities,  but  higher  dry 
densities. 

The  bituminous-stabilized  open-graded  mixtures  have  no  direct  counter- 
part in  Kersten's  series  of  tests.  Consequently,  no  comparison  of  these 
can  be  made. 

The  New  Jersey  unstabilized  open-graded  materials  in  the  moist 
condition  are  similar  to  Kersten's  crushed  granite.  In  a comparison 
with  crushed  granite  (at  120  pcf  and  U$  water  content),  the  unstabilized 
00-1  had  a conductivity  only  60$  of  that  found  by  Kersten  at  25°F,  and 
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about  80$  of  Kersten's  value  at  1*0°F.  These  figures  are  based  on  the 
Probe  test  results,  inasmuch  as  the  Probe  samples  and  Kersten's  samples 
were  both  cylindrical,  and  heat  flow  was  radial.  The  Guarded  Hot  Plate 
test  results  on  the  unstabilized  OG-1  were  h5%>  and  52%,  respectively,  of 
Kersten's  values  for  crushed  granite  at  25°F  and  Uo°F. 

A comparison  between  thermal  conductivity  values  of  the  unstabilized 
open-graded  mixtures  and  those  for  Kersten's  sandy  soils  illustrates 
that  the  equations  and  graphs  given  by  Kersten  should  not  be  used  to 
estimate  the  thermal  conductivity  of  an  open-graded  mixture.  The  Probe 
test  results  were  only  50/5  of  Kersten's  recommended  values,  both  above 
and  below  freezing;  the  Guarded  Hot  Plate  results  were  only  b0%  of 
Kersten's  values. 

A summary  of  results  from  the  open-graded  mixtures  was  given  in 
Figure  8. 


CALCULATED  FROST  DEPTHS 
Modified  Berggren  Equation 

The  modified  Berggren  equation  was  developed  by  Aldrich  and  Paynter 
(1953)  under  a contract  with  the  Arctic  Construction  and  Frost  Effects 
Laboratory  (ACFEL).  A complete  description  of  the  necessary  assumptions 
and  simplifications  made  during  its  development  are  contained  in  Aldrich 
and  Paynter  (1953).  A few  of  the  more  important  assumptions  and  limita- 
tions of  the  equation  are  discussed  below. 

Assumptions . The  mathematical  model  assumes  one-dimensional  heat 
flow  with  the  entire  soil  mass  at  its  mean  annual  temperature  prior  to 
the  initiation  of  freezing.  When  the  freezing  season  starts,  it  is 
assumed  that  the  surface  temperature  changes  as  a step  function  from  the 
mean  annual  temperature  to  the  mean  temperature  during  the  freezing 
season  and  remains  at  this  new  temperature  throughout  the  entire  freezing 
season.  The  effect  of  latent  heat  is  considered  as  a heat  sink  at  the 
moving  frost  line,  and  it  is  assumed  that  the  soil  water  freezes  at 
32.0°F. 

Limitations . The  model  on  which  the  modified  Berggren  equation  is 
based  further  assumes  an  isothermal  system  at  the  beginning  of  the 
freezing  season  and  therefore  cannot  normally  be  used  to  calculate  thaw 
depths  during  the  frost,  melting  period.  Aldrich  and  Paynter  (1953) 
state,  "Attempts  to  calculate  depth-time  curves  based  on  partial  freezing 
indices  are  likely  to  encounter  substantial  errors." 


The  modified  Berggren  equation  for  homogeneous  soils  is: 


where 

X 

K 

L 

n 

F 

X 


depth  of  freeze,  ft 

thermal  conductivity  of  the  soil,  Btu/ft  hr°F 
volumetric  latent  heat  of  fusion,  Btu/cu  ft 

conversion  factor  for  air  index  to  surface  index,  dimensionless 
air- freezing  index,  °F-days 

a coefficient  which  takes  into  consideration  the  effect  of 
temperature  changes  in  the  soil  mass 


The  X coefficient  is  a function  of  the  freezing  index,  the  mean  annual 
temperature  of  the  site,  and  the  thermal  properties  of  the  soil. 

Freezing  of  low-moisture-content  soils  in  the  lower  latitudes  is  greatly 
influenced  by  this  coefficient.  It  is  determined  by  two  factors:  the 
thermal  ratio,  a,  and  the  fusion  parameter,  u.  The  thermal  ratio  is: 


v 


v 

s 


where 

v = mean  annual  site  temperature  minus  32°F  (MAT  -32) 
o 

v = surface  freezing  index,  nF,  divided  by  the  length  of  the 
freezing  season,  t,  or  nF/t. 


The  fusion  parameter  is  equal  to  the  expression 


where 

v = same  as  previously  defined 
s 

C = average  volumetric  heat  capacity  of  the  soil, 
Btu/cu  ft°F 

L = same  as  previously  defined 


A multilayer  solution  to  the  modified  Berggren  equation  is  used  for 
nonhomogeneous , i.e.  layered,  soils,  by  determining  that  portion  of  the 
surface  freezing  index  required  to  penetrate  each  layer.  The  sum  of  the 
thicknesses  of  all  the  frozen  layers  is  the  depth  of  frost.  The  partial 
freezing  index  required  to  penetrate  the  top  layer  is  given  by 


where 


d 


1 


R 


1 


thickness  of  the  soil  layer,  ft 


hermal  resistance  of  the  layer 


sa  ft  hr  °F/Btu 


FT 


The  partial  freezing  index  required  to  penetrate  the  nth  layer  is 


where  ZR  is  the  total  thermal  resistance  above  the  nth  layer  and  equals 

ZR  = R,  + R,,  + R_  . . . + R , 

123  n-1 

When  the  summation  of  the  partial  indexes,  F + + F . ..+  F , is  equal 

to  the  surface  freezing  index,  the  frost  depth  has  been  determined. 


Aitken  and  Berg  (1968)  presented  a computer  program  solving  the 
modified  Berggren  equation  for  multiple-layered  soil  and  pavement  pro- 
files. A revised  version  of  that  computer  program  was  used  in  this 
study;  this  program  is  listed  in  Appendix  D.  Aitken  and  Berg  discuss, 
in  greater  detail,  the  various  parameters  that  are  used  in  the  computer 
program  to  solve  the  modified  Berggren  equation. 


Comoutat ions 


To  establish  initial  and  upper  boundary  conditions  for  the  modified 
Berggren  equation  solutions , monthly  and  annual  weather  records  from 
weather  stations  in  New  Jersey  were  surveyed.  Three  sites  with  long 
term  records  were  chosen  for  the  calculations.  Trenton  and  Atlantic 
City  were  chosen  for  southern  sites  and  Newton  was  chosen  as  representa- 
tive of  a northern  site.  Data  used  in  the  modified  Berggren  equation 
solutions  are  summarized  in  Table  8. 


Table  8.  Initial  and  boundary  conditions  for 
modified  Berggren  equation  solutions. 


Mean  Freezing  Season 

Atlantic 

Newton 

Trenton 

City 

Air  Freezing  Index,  °F-days 

1+80 

130 

50 

Length  of  Season,  days 

95 

1*0 

10 

Mean  Annual  Air  Temperature,  °F 

L9.6 

53.8 

5l*  - 1* 

M-Factor  (Portland  cement) 

O.M* 

0.1+1* 

0.1*1* 

N-Factor  (Asphaltic  concrete) 

0.50 

0.50 

0.50 

Design 

Freezing  Season 

Atlantic 

Newton 

Trenton 

City 

Air  Freezing  Index,  °F-days 

900 

320 

250 

Length  of  Season,  days 

110 

83 

82 

Mean  Annual  Air  Temperature,  °F 

L9.6 

53.8 

5I+.I+ 

N-Factor  ( Portland  cement ) 

0.1*1* 

0.1*1* 

0.1*1* 

N-Factor  (Asphaltic  concrete) 

0.50 

0.50 

0.50 
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Two  sets  of  freezing  season  data  were  used  in  the  solutions,  one  an 
average  or  mean  freezing  season  and  the  other  a design  freezing  season. 

The  design  freezing  season  is  chosen  as  the  average  of  the  three  coldest 
winters  in  the  1946-1976  period. 

Table  9 contains  gradation  specifications  for  base  and  subbase 
materials  used  in  New  Jersey  highway  pavements , Table  10  contains  thick- 
nesses of  various  layers  in  typical  New  Jersey  pavement  sections  without 
a drainage  layer,  and  Table  11  contains  proposed  thicknesses  of  various 
layers  for  New  Jersey  pavement  sections  incorporating  either  stabilized 
or  unstabilized  drainage  layers.  Both  of  these  sets  of  information  were 
furnished  by  engineers  from  the  New  Jersey  Department  of  Transportation. 
Thermal  properties  for  pavement  sections  without  drainage  layers  are  given 
on  Table  12  and  thermal  properties  for  pavement  cross  sections  incorporating 
drainage  layers  on  Table  13. 


Table  9.  Gradation  specifications  for  base  and  subbase 
materials  in  New  Jersey 


Designation 

1A 

IB 

1C 

2B 

5A 

Sieve  Size 

Percentage  by 

Weight  F 

assing 

4" 

100 

100 

2-1/2" 

100 

2" 

70-100 

100 

100 

3/4" 

50-95 

65-100 

60-100 

70-100 

50-90 

No.  4 

30-60 

40-75 

30-100 

30-80 

25-60 

No.  50 

5-25 

5-30 

5-35 

10-35 

5-25 

No.  200 

0-7 

0-7 

0-5 

5-12 

3-12 

NOTES:  1A  is  usually  used  as  a base  course  beneath  PCC  pavement. 


IB  may  be  used  as  a base  and/cr  subbase  course  beneath  a PCC 
pavement . 

1C  is  generally  used  as  a subbase  beneath  1A  or  5A  base  courses . 

2B  is  a gravel  base  course  commonly  used  beneath  shoulders  and 
secondary  roads . 

5A  is  a quarry -processed  material  and  is  used  beneath  bituminous 
stabilized  base  course,  in  bituminous  concrete  pavements. 


w - 

^ 

Table  10.  Typical  New  Jersey  pavement  sections 
without  a drainage  layer. 

f 

n 

. 

Section  #1  Southern  New  Jersey: 

Section  H 3 Northern  New  Jersey: 

Section  #5  Southern  New  Jersey: 

Section  #6  Southern  New  Jersey 

Section  #7  Northern  New  Jersey 

PCC  - Portland  Cement  Concrete 
FABC  - Fine  Aggregate  3ituminous  Course 
3SBC  - Bituminous  Stabilized  3ase  Course 

•Values  for  1C  material  were  used  in  the  thermal  analysis. 
See  Tables  6 and  9 for  properties  of  the  granular  materials. 


9"  PCC 

12"  Base  (IB) 
Subgrade* 

9"  PCC 
6"  Base  (1A) 

6"  Subbase  (1C) 
Subgrade* 


3"  FABC 
6"  BSBC 
5"  3ase  (2B) 

5"  Subbase  (13) 


3"  ~ABC 
6"  3SBC 
6"  Base  (5A) 

11"  Subbase  (IB) 


3"  FABC 
6"  BSBC 
6"  Base  (5A) 

12"  Subbase  (1C) 
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Table  11. 


Tentative  New  Jersey  pavement  sections 
incorporating  a drainage  layer . 


Portland  Cement  Concrete  Pavements 


Section  #2  Southern  New  Jersey 
2A  OGL  is  stabilized 
2B  OGL  is  unstabilized 


6"  drainage  layer  (OGL) 

6"  Subbase  (lB  or  2B)  with  the 
top  2"  lime  fly  ash  or  cement 
stabilized. 


Section  Northern  New  Jersey 
UA  OGL  is  stabilized 
UB  OGL  is  unstabilized 


9"  PCC 

6"  Drainage  layer  (OGL) 
6"  Subbase  (lA)  with  2" 
Stabilization 


Bituminous  Concrete  Pavements 


3"  FABC 

6"  BSBC  course 
5"  Drainage  layer  (OGL) 

2"  Subbase  (IB)  Stabilized 
w/cement  or  lime-fly  ash 
3"  Subbase  (IB) 


3"  FABC 
6"  BSBC 

6"  Drainage  layer  (OGL) 

2"  Subbase  (1C)  stabilized 
w/cement  or  lime-fly  ash 
10"  Subbase  (1C) 

GL  - Open  graded  layer.  This  layer  may  be  stabilized  with  an  asphalt 
cement  in  either  type  of  pavement  or  used  unstabilized  under  Port- 
land cement  concrete  pavements. 
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Section  #8  Southern  New  Jersey 


Section  #9  Northern  New  Jersey 


k 


Table  12.  Thermal  properties  for  pavement  sections 
without  a drainage  layer 


SECTION  1 


PCC 

lUo 

0 

0.75 

0.57 

_ 

30.0 

0 

Base  ( IB) 

129 

6. 7(9- 7) 

1 .00 

1.15 

1.U5 

28.  Ul 

31.31 

12U5 

1802 

Subgrade  ( 1C ) 

118 

lU.U 

5.00 

- 

1. 08 

- 

32.80 

- 

2UU7 

SECTION  3 

PCC 

lUo 

0 

0.75 

0.57 

• 

30.0 

_ 

• 

Base  (1A) 

129 

6. 7(9- 7) 

0.50 

1.15 

1.1*5 

28.  Ul 

31.31 

12U5 

1082 

Subbase  (1C) 

117 

6.3(12.0) 

0.50 

1.17 

1.81 

25. U2 

30. U2 

1061 

2022 

Subgrade  (lC) 

118 

lU.U 

10.00 

- 

1.08 

- 

32.80 

- 

2UU7 

SECTION  5 

FAEC 

138 

0 

0.25 

0.5U 

— 

28.0 

— 

0 

- 

BS3C 

138 

0 

0.  50 

0.53 

- 

28.0 

- 

0 

_ 

Base  (23) 

129 

6.7(9. 7) 

0.U2 

1.15 

1.1*5 

28.  Ul 

31.31 

12U5 

1802 

Subbase  (IB) 

129 

6. 7(9. 7) 

0.U2 

1.15 

1.U5 

28.  Ul 

30. U2 

12U5 

1802 

Subgrade  ( 1C ) 

118 

lU.U 

5.00 

- 

1.06 

- 

32.80 

- 

2UU7 

SECTION  6 

FABC 

138 

0 

0.25 

0.5U 

28.0 

— 

0 

- 

BSBC 

138 

0 

0.50 

0.53 

- 

28.0 

- 

0 

- 

Base  (5A) 

128 

8.1(11.7) 

0.  50 

0.99 

1.05 

29- 5U 

37. 5U 

1U93 

2157 

Subbase  (13) 

129 

6. 7(9. 7) 

0.92 

1.15 

1.U5 

28.  Ul 

31.31 

12U5 

1802 

Subgrade  ( 1C ) 

118 

lU.U 

15.00 

- 

1.08 

- 

32.80 

- 

2UU7 

SECTION  T 

FA3C 

138 

0 

0.25 

0.5U 

* 

28.0 

0 

BSBC 

138 

0 

0.50 

0.53 

- 

28.0 

- 

0 

- 

Base  (5A) 

128 

8.1(11.7) 

0.50 

0.99 

1.05 

29.  5U 

32.99 

1U93 

2157 

Subbase  ( 1C ) 

117 

6.3(12.0) 

1.00 

1.17 

1.8l 

25. U2 

30. U2 

1061 

2022 

Note:  Moisture  contents  in  parentheses  are  the  "wet"  values. 


Table  13.  Thermal  properties  for  pavement 
sections  with  a drainage  layer 


Material 

Yd 

W 

d 

k 

C 

L 

lb/ftJ 

% 

ft 

Btu 

fthr°F 

Btu 
ft  ioF 

Btu 

ft7 

SECTION  2A 

PCC 

lUo 

0 

0.75 

0.57 

30.0 

0 

OGL  (stabilized) 

120 

7.0 

0.5 

0.51 

26.6 

1201* 

Subbase  (lB) 

129 

6.7 

0.5 

1.15 

28.1* 

12l*5 

Subgrade  ( 1C ) 

118 

Ih.h 

5.0 

1.08 

32.8 

21*1*7 

SECTION  2B 

PCC 

lUO 

0 

0.75 

0.57 

30.0 

0 

OGL  (unstabilized) 

121 

h.5 

0.50 

0.1*0 

25.9 

781* 

Subbase  ( IB ) 

129 

6.7 

0.50 

1.15 

28. 1* 

121*5 

Subgrade  ( 1C ) 

118 

Ih.h 

5.00 

1.08 

32.8 

2U1*7 

SECTION  1*A 

PCC 

lUO 

0 

0.75 

0.57 

30.0 

0 

OGL  (stabilized) 

120 

7.0 

0.5 

0.51 

2 6.6 

1201* 

Subbase  (lA) 

129 

6.7 

0.5 

1.15 

28.1* 

121*5 

Subgrade  (1C) 

118 

1U.U 

5.0 

1.08 

32.8 

2UL7 

SECTION  1*3 

PCC 

lUO 

0 

0.75 

0.57 

30.0 

0 

OGL  (unstabilized) 

121 

h.5 

0. 50 

0.1*0 

25.9 

781* 

Subbase  (lA) 

129 

6.7 

0.50 

1.15 

28.1* 

12l*5 

Subgrade  ( 1C ) 

118 

ll*.l* 

5.00 

1.08 

32.8 

21*1*7 

SECTION  8 

FABC 

138 

0 

0.25 

0.51* 

28.0 

0 

BSBC 

138 

0 

0.50 

0.51 

28.0 

0 

OGL  (stabilized) 

120 

7.0 

0.1*2 

0.51 

2 6.6 

1201+ 

Subbase  (13) 

129 

6.7 

0.1*2 

1.15 

28.1* 

121*5 

Subgrade  ( 1C ) 

118 

Ih.h 

5.0 

1.08 

32.8 

21*1*7 

SECTION  9 

FABC 

138 

0 

0.25 

0.5l* 

28.0 

0 

BSBC 

138 

0 

0.50 

0.5I* 

28.0 

0 

OGL  (stabilized) 

120 

7.0 

0.50 

0.51 

2 6.6 

1201* 

Subgrade  ( 1C ) 

118 

lU.l* 

5.0 

1.08 

32.8 

21*1*7 
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Values  of  thermal  properties  snown  in  the  tables  are  average  values 
for  frozen  and  unfrozen  conditions.  The  volumetric  heat  capacity,  which 
is  the  quantity  of  heat  required  to  change  the  temperature  of  a unit 
volume  by  one  degree,  is  determined  from 

C = Yd  [c  + 0.75  (w/100 ) J 


where 

C = volumetric  heat  capacity,  Btu/cu  ft°F 
Yd  = dry  unit  weight  of  soil,  lb/cu  ft 
c = specific  heat  of  dry  soil,  Btu/lb  °F 
w = moisture  content  of  soil,  % dry  weight. 

The  volumetric  latent  heat  of  fusion  is  the  quantity  of  heat  re- 
quired to  freeze  the  water  in  a unit  volume  of  soil  without  a change  in 
temperature.  It  is  determined  from 

L = lUU  y (w/100) 
d 

where  the  symbols  are  as  defined  above. 

The  thermal  conductivity  values  used  in  most  of  the  computations 
were  determined  from: 

Ka  - K*Kf)/2 


where 

K = average  thermal  conductivity,  Btu/ft  hr°F 

3. 

= thermal,  conductivity  at  about  U0°F 

Kj,  = thermal  conductivity  at  about  22°F 

The  values  of  K and  were  generally  those  obtained  from  the 
laboratory  measurements  made  during  this  study.  Thermal  conductivity 
values  of  the  Portland  cement  and  asphaltic  concrete  specimen  furnished 
to  CRREL  by  the  NJDOT  were  lower  than  the  values  normally  used.  There- 
fore several  frost  depth  calculations  were  made  using  the  larger  values. 
The  outcome  of  the  computations  is  discussed  in  the  section  "Results." 

Information  from  Tables  8,  10,  11,  12  and  13  was  used  in  the  com- 
puter program  discussed  previously  to  calculate  maximum  frost  penetra- 
tion depths  for  various  conditions.  For  each  pavement  profile  shown  in 
Table  10  solutions  were  run  using  the  mean  freezing  index  and  the  thermal 
properties  in  the  dry  state  to  compute  a maximum  frost  penetration  depth 
for  the  mean  year  and  again  in  the  wet  state  to  comrute  a minimum  frost 
penetration  depth  for  the  mean  year.  Two  similar  sets  of  calculations 
were  made  for  each  of  the  design  freezing  seasons.  Therefore,  four 
frost  penetration  depths  were  computed  for  each  of  the  pavement  sections 
listed  in  Table  10.  Calculations  were  made  using  only  the  "dry"  soil 
properties  for  profiles  shown  in  Table  11. 
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Results 


Table  lU  contains  the  results  of  the  computations.  The  first  part 
of  the  table  contains  data  from  pavement  profiles  without  drainage 
layers  and  the  second  part  contains  data  for  pavement  profiles  incorporating 
drainage  layers.  Frost  penetration  in  Section  1 progresses  only  slightly 
beneath  the  Portland  cement  concrete  pavement.  In  Section  3 the  frost 
penetrates  slightly  into  the  base  material,  and  in  Section  5,  the  frost 
barely  penetrates  the  pavement  and  the  bituminous  stabilized  base  course. 

In  Section  6 the  same  phenomenon  occurs  and  in  Section  7 frost  penetra- 
tes into  the  base  course  beneath  the  asphalt  pavement  and  the  bituminous 
stabilized  base  course.  Thus,  for  pavement  profiles  without  drainage 
layers  the  calculated  maximum  frost  depth  in  the  southern  part  of  the 
state  barely  penetrates  the  pavement  section  during  either  a mean  or  a 
design  season  or  with  wet  or  dry  subsoils.  In  the  northern  part  of  Hew 
Jersey,  frost  may  penetrate  into  the  subbase  material. 

The  last  six  calculations  in  the  first  part  of  Table  lH  were  made 
to  illustrate  the  influence  of  higher  thermal  conductivities  of  the 
pavement  layer.  For  the  southern  sites,  increasing  the  thermal  conduc- 
tivity of  the  pavement  did  not  significantly  alter  the  computed  frost 
depths.  At  the  northern  site,  the  conputed  frost  depth  increased  by 
about  25$  when  the  thermal  conductivities  of  the  pavement  layers  were 
increased  to  the  values  shown  in  the  footnote  to  Table  lk. 

Data  in  the  lower  part  of  Table  lU  indicate  that  frost  penetration 
into  pavement  profiles  with  open-graded  layers  is  equal  to  or  slightly 
less  than  beneath  similar  pavement  profiles  without  drainage  layers. 

The  difference  is  due  to  the  lower  thermal  conductivity  of  the  open-graded 
material.  Whether  the  open-graded  material  is  stabilized  or  unstabilized 
has  little  or  no  influence  on  the  depth  of  frost  penetration  into  the 
pavement . 

These  computations  suggest  that  frost  penetration  beneath  a New 
Jersey  highway  containing  a drainage  layer  will  be  equal  to  or  slightly 
less  than  the  frost  penetration  beneath  a "typical"  highway  currently 
constructed  there. 

The  n-factors  shown  on  Table  8 were  used  to  compute  the  maximum  frost 
depths  shown  on  Table  lU.  Had  the  average  n-f actor  for  Portland  cement 
and  asphaltic  concrete  pavements  shown  on  Table  1 been  used  in  the  cal- 
culations, maximum  frost  depths  on  Table  lU  would  increase  by  8$  for 
Portland  cement  concrete  pavements  and  13$  for  asphaltic  concrete  pave- 
ments. Thus  the  maximum  difference  in  the  calculated  values  would  be  an 
increase  of  0.3  ft  for  Section  7 **  in  Table  lU. 

Frost  penetration  depths  calculated  in  this  report  have  used  a design 
freezing  index  which  was  determined  by  the  average  of  the  three  coldest 
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Table  lU.  Maximum  seasonal  frost  penetration  determined 
with  the  modified  Berggren  equation. 


Profiles  without  drainage  layers : 

Maximum  Seasonal 

Frost  Penetration  Thickness  of 


Section 

Conditions 

Dry  Soil 

Wet  Soil 

pavement  section 

1 

Mean  Season  (Atlantic  City) 

0.8  ft 

0.8  ft 

1.75  ft 

1 

Mean  Season  (Trenton) 

0.8 

0.8 

1.75 

1 

Design  Season  (Atlantic  City) 

0.8 

0.8 

1.75 

1 

Design  Season  (Trenton) 

0.8 

0.8 

1.75 

3 

Mean  Season  (Newton) 

1.0 

0.9 

1.75 

3 

Design  Season  (Newton) 

1.6 

1.6 

1.75 

5 

Mean  Season  (Atlantic  City) 

0.8 

0.8 

1.58 

5 

Mean  Season  (Trenton) 

0.8 

0.8 

1.58 

5 

Design  Season  (Atlantic  City) 

0.8 

0.8 

1.58 

5 

Design  Season  (Trenton) 

0.8 

0.8 

1.58 

6 

Mean  Season  (Atlantic  City) 

0.8 

0.8 

2.17 

6 

Mean  Season  (Trenton) 

0.8 

0.8 

2.17 

6 

Design  Season  (Atlantic  City) 

0.8 

0.8 

2.17 

6 

Design  Season  (Trenton) 

0.8 

0.8 

2.17 

7 

Mean  Season  (Newton) 

1.0 

1.0 

2.25 

7 

Design  Season  (Newton) 

1.7 

1.5 

2.25 

1* 

Design  Season  (Atlantic  City) 

0.8 

- 

1.75 

1* 

Design  Season  (Trenton) 

0.8 

- 

1.75 

3* 

Design  Season  (Newton) 

2.0 

- 

1.58 

5** 

Design  Season  (Atlantic  City) 

0.8 

- 

1.58 

5## 

Design  Season  (Trenton) 

0.8 

- 

1.58 

y## 

Design  Season  (Newton) 

2.1 

- 

2.25 

Profiles 

with  drainage  layers : 

Thickness  of 

Section 

Conditions 

Dry  Soil 

pavement  section 

2A 

Mean  Season  (Atlantic  City) 

0.8  ft 

1.75  ft 

2A 

Mean  Season  (Trenton) 

0.8 

1.75 

2A 

Design  Season  (Atlantic  City) 

0.8 

1.75 

2A 

Design  Season  (Trenton) 

0.8 

1.75 

23 

Mean  Season  (Atlantic  City) 

0.8 

1.75 

2B 

Mean  Season  (Trenton) 

0.8 

1.75 

2B 

Design  Season  (Atlantic  City) 

0.8 

1.75 

2B 

Design  Season  (Trenton) 

0.8 

1.75 

UA 

Mean  Season  (Newton) 

1.0 

1.75 

UA 

Design  Season  (Newton) 

l.U 

1.75 

UB 

Mean  Season  (Newton) 

0.9 

1.75 

US 

Design  Season  (Newton) 

1.3 

1.75 
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Table  ll+ 

(cont’ d) 

Thickness 

Section 

Conditions 

Dry  Soil 

pavement  section 

8 

Mean  Season  (Atlantic  City) 

0.8 

1.58 

8 

Mean  Season  (Trenton) 

0.8 

1.58 

8 

Design  Season  (Atlantic  City) 

0.8 

1.58 

8 

Design  Season  (Trenton) 

0.8 

1.58 

9 

Mean  Season  (Newton) 

1.0 

1.1+2 

9 

Design  Season  (Newton) 

1.1+ 

1.1+2 

* Thermal  conductivity  of  PCC  increased  to  1.0  Btu/ft  hr  °F,  other  data  same  as 
in  Table  12. 

**  Thermal  conductivities  of  FABC  and  BSBC  increased  to  0.86  Btu/ft  hr  °F,  other 
data  same  as  in  Table  12. 

winters  in  the  last  30  years  of  record.  Figure  14,  from  Department  of 
the  Army  (1966),  illustrates  the  relationship  between  mean  air  freezing 
indexes  and  indexes  determined  for  colder  years  for  30  consecutive 
years.  Another  method  of  developing  design  air  freezing  index  values  is 
to  use  the  relationships  developed  by  Berube  (1967),  Figure  15*  His 
method  is  statistically  based  and  the  designer  can  choose  the  design 
life  and  the  level  of  probability  that  a computed  freezing  index  will  be 
exceeded  during  the  design  period.  The  computed  freezing  index  can  then 
be  used  to  estimate  the  frost  depth.  An  example  of  the  relationships 
that  can  be  obtained  from  this  procedure  is  shown  in  Figure  16.  Conditions 
similar  to  those  from  Newton  where  the  mean  freezing  index  equals  l+80°F-days 
(Table  8),  with  a pavement  cross  section  similar  to  Section  7 were  used 
in  the  computations  for  the  figure.  A design  engineer  may  wish  to  study 
the  effect  of  differing  design  life  on  the  maximum  frost  penetration 
depth  beneath  the  pavement  surface.  He  may  also  wish  to  evaluate  the 
effect  on  frost  penetration  of  different  probabilities  of  a design  index 
being  exceeded.  Results  from  a sample  study  are  shown  in  Table  15.  These 
indicate  that  for  a given  design  life,  the  frost  depth  increases  as  the 
probability  that  the  design  freezing  index  will  be  exceeded  decreases. 

Note  that  the  probability  the  900  °F-day  design  index  previously  used  in 
the  calculations  will  be  exceeded  is  85$  for  a 30-year  design  life  but 
only  50$  for  a 10-year  design  life.  The  frost  depth  varied  by  more  than 
h0%  (2. U-I.7/I.7)  depending  upon  the  design  life  and  the  probability  that 
the  design  freezing  index  will  be  exceeded  at  least  once  during  this 
design  life.  Using  Figure  15  the  design  engineer  can  determine  a design 
frost  depth  based  on  the  design  life  and  importance  of  the  road,  i.e.  the 
more  important  the  road,  the  lower  the  allowed  probability  that  the  design 
freezing  index  can  be  exceeded. 

OTHER  CONSIDERATIONS 

The  coldest  winter,  which  will  normally  cause  the  deepest  frost 
penetration,  may  not  cause  the  most  extensive  damage  to  pavements.  For 
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Fisrure  l6.  Influence  of  freezing  index  on  frost  pene- 
tration beneath  Section  7 . 
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Figure  15.  Statistically  based  design  air  treezing  index  values 
( Berube , 1967 ) . 


Table  15. 

Thaw  depths  from  a 

statistically  based 

design  method 

Design 

Probability  of 

Design 

Frost 

Life 

Being  Exceeded 

Index 

Depth 

Years 

Of 

JO 

°F-days 

ft 

10 

5 

1120 

2.1 

10 

10 

1080 

2.0 

10 

20 

1020 

1.9 

10 

50 

900 

1.7 

30 

5 

1260 

2.4 

30 

10 

1180 

2.2 

30 

20 

1130 

2.1 

30 

50 

1020 

1.9 

30 

«5 

900 

1.7 

50 

5 

1280 

2.4 

50 

10 

1230 

2.3 

50 

20 

1170 

2.2 

50 

50 

1070 

2.0 

example,  Cosaboom  (1977)  noted  that  the  1975-1976  winter,  although  rela- 
tively mild,  caused  several  freeze-thaw  cycles  in  the  upper  pavement 
layers,  resulting  in  extensive  damage  to  New  Jersey  pavements.  In  con- 
trast the  1976-1977  winter  was  extremely  cold  with  very  few  freeze-thaw 
cycles  in  the  pavement  layers.  Relatively  little  pavement  damage  was 
experienced  during  the  1976-1977  winter. 

Since  quantitative  data  on  moisture  contents,  precipitation  and 
temperatures  within  the  pavement  systems  are  not  available,  it  is 
impossible  to  state  with  certainty  the  exact  cause  of  this  difference  in 
behavior.  From  the  qualitative  data  which  are  available,  and  from 
results  of  field  and  laboratory  tests  conducted  by  others,  for  example 
Jessberger  and  Carbee  (1970)  and  Fredrickson  (1963),  on  the  reduction  of 
strength  of  soils  upon  thawing,  the  major  problem  is  caused  by  excessive 
water  in  the  previously  frozen  material.  Frequent  freeze-thaw  cycles 
may  permit  more  opportunity  for  moisture  to  enter  the  subsurface  materials 
through  cracks  in  the  pavement.  If  vertical  and  lateral  drainage  do 
not  permit  sufficient  water  removal  between  successive  freeze-thaw 
cycles,  expansion  of  excess  water  upon  freezing  can  damage  pavements  at 
cracks  and  joints,  and  extensive  weakening  at  the  cracks  and  joints  may 
cause  additional  breakage  upon  thawing.  Incorporation  of  a drainage 
layer  beneath  the  pavement  should  reduce  the  severity  of  this  problem. 

The  drainage  layer  must  be  part  of  a reliable  subsurface  drainage  system 
wherein  all  components  function  satisfactorily  under  frozen  and  unfrozen 
conditions.  Cedegren  et  al.  (1972),  the  Department  of  Transportation 
(1972)  and  Cedegren  (1974)  provide  guidelines  for  designing  and  in- 
stalling drainage  layers  for  pavement  systems . 
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Since  the  thermal  conductivities  of  the  open-graded  drainage  layers 
are  generally  somewhat  lower  than  the  thermal  conductivities  of  the  base 
and  subbase  course  materials,  a question  has  been  raised  as  to  whether 
this  will  cause  the  pavement  surface  over  an  open-graded  drainage  layer 
to  become  icy  before  a pavement  without  the  drainage  layer  does  (Cosaboom, 
1977).  This  problem  has  occurred  when  insulation  layers  are  incorporated 
in  pavement  cross-sections  and  is  similar  to  the  behavior  of  bridge 
decks  which  sometimes  become  ice-covered  when  the  adjacent  roadway  is  not. 
The  difference  in  thermal  conductivity  between  the  open-graded  layer 
(0.5^  Btu/ft  hr  °F)  and  the  base  course  and  subbase  course  materials 
(about  1.1  3tu/ft  hr  °F)  is  considerably  less  than  the  difference  between 
a thermal  insulating  layer  (0.02  Btu/ft  hr  °F)  and  the  base  course  and 
subbase  course  materials.  Therefore  we  do  not  anticipate  any  significant 
difference  in  surface  conditions  between  pavements  with  and  without  an 
open-graded  drainage  layer.  Other  environmental  conditions,  i.e.  shaded 
or  open,  south  side  or  north  side  of  a hill,  etc.,  will  probably  cause 
more  differences  in  surface  conditions  than  the  use  of  an  open-graded 
drainage  layer.  We  suggest,  however,  that  when  segments  of  highways 
containing  open-graded  drainage  layers  are  constructed  in  New  Jersey, 
careful  observations  be  made  to  determine  whether  differential  icing 
conditions  occur  over  these  pavements. 


FUTURE  WORK 

At  the  suggestion  of  CRREL  engineers,  New  Jersey  Department  of 
Transportation  engineers  and  technicians  installed  frost  tubes  for 
monitoring  frost  penetration  depths  and  access  tubes  for  use  in  monitoring 
subsurface  moisture  contents  with  a nuclear  probe.  During  the  1976-1977 
winter,  the  New  Jersey  Department  of  Transportation  measured  frost 
penetration  depths  and  moisture  contents  at  several  sites.  Air  and 
surface  temperatures  will  again  be  measured  at  selected  sites  in  New 
Jersey  and  n-factors  and  freezing  indexes  will  be  calculated  from  the 
temperature  data.  The  New  Jersey  Department  of  Transportation  has 
furnished  these  data  to  CRREL  and  we  are  comparing  measured  frost  pene- 
tration depths  with  depths  calculated  using  the  modified  Berggren 
equation.  These  data  will  be  used  to  evaluate  our  confidence  in  the 
calculated  frost  penetration  depths  and  to  validate  our  tentative  con- 
clusion that  the  use  of  drainage  layers  should  have  no  detrimental 
effects  on  frost  action  in  New  Jersey  pavement  systems. 
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APPENDIX  A 


I 


Material  Characteristics  Furnished  by  the  NJDOT 


Table  Al.  Specifications  for  Bituminous  Concrete  Mixtures 


k 


Gradation 

Mix  No. 

Base 

1 

Binder 

2 

MA-BC 

h 

FA-BC 

5 

Sieve 

Size 

Grading  of  total  aggregate  (coarse  plus  fine,  plus 
filler  if  required).  Amounts  finer  than  each  lab- 
oratory sieve  (square  opening),  weight  percent. 

2" 

100 

— 

— 

— 

1-1/2" 

90-100 

100 

— 

— 

1" 

80-100 

90-100 

— 

— 

3 A" 

— 

— 

100 

1/2" 

50-85 

60-80 

90-100 

100 

3/8" 

— 

— 

— 

80-100 

No.  1* 

25-60 

25-60 

U0-70 

55-75 

No.  k 

20-50 

15  A 5 

25-55 

30-60 

No.  1 6 

— 

— 

— 

— 

No.  30 

— 

— 

— 

— 

No.  50 

8-30 

3-18 

10-25 

10-30 

No.  100 

— 

— 

— 

— 

No.  200 

U-12 

1-7 

H-10 

1+-10 

NOTE:  Material  Passing  The  No.  200  Sieve  May  Consist  of 
Fine  Particles  of  the  Aggregate  or  Mineral  Filler,  or  Both.  In 
Accordance  with  the  Requirements  of  Current  A.A.S.H.O.  Designation 
T 90. 


Table  A2.  Bituminous  Mix  #1. 


Materials : 

#57  Stone  - Houdaille  Construction  Materials,  Bound  Brook,  N.J. 
(Diabase ) 

#8  Stone  - Houdaille  Construction  Materials,  Bound  Brook,  N.J. 
(Diabase ) 

Sand  - Houdaille  Construction  Materials,  Kenville,  N.J. 

Mineral  Filler  - Limestone  Products,  LimeCrest , N.J. 

Asphalt  Cement  - ARCO  AC-20 


II. 

Aggregate 

Gradat ion 

#57 

#8 

Sand 

Sieve 

Specs . 

Actual 

Size 

{%  uassing) 

(<* 

\ fa 

nassing ) 

Suecs . 

Actual 

Specs.  Actual 

2” 

100 

100 

100 

100 

100 

1-1/2" 

100 

100 

100 

100 

1" 

95-100 

99.7 

100 

100 

100 

3 A" 

88.0 

100 

100 

100 

1/2" 

25-60 

1+7-5 

100 

100 

100 

3/8" 

29.3 

85-100 

96.8 

100 

#1+ 

0-10 

7.0 

10-30 

23.1 

98.1 

#8 

0-5 

2.2 

0-10 

3.1 

88.1 

#50 

2.1 

22.7 

#200 

1.5 

III.  Cylinder  Properties 

A. 

Voids  - 6 

(in 

each ) 

3. 

Dimensions  - 

1) 

6"  (Diam. 

) x 8.30" 

2) 

6"  (Diam. 

) x 8.1+2" 

C. 

Weights  of  materials  used 

( each ) 

1)  #57 

_ 

3LL0 

gms . 

2)  #8 

— 

2150 

gUlS  • 

3 ) Sand  — 

— 

2150 

gms . 

1+ ) Mineral 

Filler 1+lU 

gms  • 

5 ) Asphalt 

Cement  — 1+53 

gms . 

Table  A2  (cont'd) 


IV. 


Cylinders  Design  Specifications 
Bituminous  Mix  #1  - Base  Course 


Sieve  Size 

Specifications  ($  passing) 

Actual  (,%  passing) 

2" 

100 

100 

1-1/2" 

90-100 

100 

1" 

80-100 

99.8 

3A" 

95.2 

1/2" 

50-85 

79.0 

3/8" 

70.7 

#U 

25-60 

39.2 

#8 

20-50 

28.8 

#50 

8-30 

11.3 

#200 

U-12 

5.2 

Asphalt  Cement , Weight  Percentage  of  Total  Mixtures : 


Spec. 3. 5-8.0 

Actual 5.0 


Table  A3.  Bituminous  Mix  #2  - Binder  Course 


Materials 

#57  Stone  - Houdaille  Construction  Materials,  Bound  Brook,  N.J. 
(Diabase ) 

# 8 Stone  - Houdaille  Construction  Materials,  Bound  Brook,  N.J. 
(Diabase) 

Sand  - Houdaille  Construction  Materials,  Kenville,  N.J. 

Asphalt  Cement  - Arco  AC-20 


II.  Aggregate  Cradation 


#57 


#8 


Sand 


Sieve 

Specs . 

Actual 

Soecs . 

Actual  Specs . 

Actual 

Size 

(%  passing) 

[f  passing) 

( % passing) 

(%  passing) 

2" 

100 

100 

100 

100 

100 

1-1/2" 

100 

100 

100 

100 

100 

1" 

95-100 

97.7 

100 

100 

100 

3/U" 

100 

100 

100 

1/2" 

25-60 

A.  5 

100 

100 

100 

3/8" 

85-100 

98.3 

100 

#1* 

0-10 

0.8 

10-30 

35.6 

99.1 

#8 

0-5 

0.7 

0-10 

6.7 

38.1 

#50 

1.2 

16.5 

#200 

1.2 

1.9 

!*5 


'able  A3  (cont'd) 


III.  Physical  Properties 

Once  the  pre-measured  material  (aggregate,  sand,  and  asphalt  cement) 
has  been  combined  and  compacted  into  a 20"  x 20"  x 3"  form,  one 
should  have  a representative  #2  mixture  with  67  voids 

Weights  of  materials  used 

1)  #57 20, 577. U gms 

2)  #8 13,718,3  gms 

3)  Sand 11,1+31.9  gms 

U)  Asphalt  Cement....  2,1+06.7  gms 


TV.  Bituminous  Mix  #2  (Binder  Course) 


Sieve  Size 

Specs.  (7  passing) 

Actual  (7  passing) 

2" 

1-1/2" 

100 

100 

1" 

90-100 

99-2 

3/1+ " 

1/2" 

60-80 

70.1 

3/8" 

#1+ 

25-60 

33.0 

#8 

15-1+5 

16.9 

#50 

3-13 

3.1 

#200 

1-7 

1.0 

Asphalt  Cement, 

Weight  Percentage  of  Total 

Mixture : 

Specs 

1+.0  - 3.57 

Actual. . . . 

5.07 

Table  AU.  Bituminous 

Mix  #5. 

T.  Materials: 


#8  Stone  - Houdaille  Construction  Materials,  Bound  Brook,  N..T. 
(Diabase ) 

Sand  - Houdaille  Construction  Materials,  Kenville,  N.J. 

Mineral  Filler  - Limestone  Products,  Limecrest,  N.J. 

Asphalt  Cement  - ARCO  AC-20 


1+6 


Table  A4  ( cont ' d ) 


II. 


Aggregate-  Gradation 


Sieve  Size 

Soecs  (%  pass) 

2" 

100 

1-1/2" 

100 

1" 

100 

3/4" 

100 

1/2" 

100 

3/8" 

85-100 

#4 

10-30 

#6 

0-10 

#50 

#200 

Actual  Specs . Actual 


100 

100 

96.8 

100 

23.1 

98.1 

3.1 

88.1 

2.1 

22.7 

1.5 

III.  Cylinder  Properties 


A. 

Voids  6%  (each) 

3. 

Dimensions 

1)  6" 

1 ( diam. ) x 

8.26" 

2)  6” 

' (diam. ) x 

8.40" 

C. 

Weights  of  : 

materials 

used 

1) 

#8-  — 

37  64 

gms 

2) 

Sand  — 

1207 

gms 

3) 

Mineral 

Filler  — 

421  gms 

4) 

Asphalt 

Cement  — 

443  gms 

IV.  Cylinder's  Design  Specifications 


3ituminous  Mix  #5  - Surface  Course 


Table  A5.  Specific  Gravity  of  Bituminous  Cylinder  Materials 


# 57 2.91 

#8  Stone 2.93 

Sand 2.63 


Mineral  Filler- 
Asphalt  Cement 


2.75 

1.01 


Table  A 6.  Portland  Cement  Concrete  Cylinders 

I.  Materials: 

Portland  Concrete  Cement 
Concrete  Sand  (Fine)  Aggregate 
Coarse  Aggregate 


II.  Aggregate  Gradations  {%  passing  by  weight) 


Fine 

Aggregate 

Coarse 

Aggregate 

Sieve  Size 

Soecs . 

Actual 

Soecs . 

Actual 

1-1/2" 

100 

100 

1" 

95-100 

100 

1/2" 

25-60 

32 

3/8" 

100 

100 

#U 

95-100 

99 

0-10 

h 

# 8 

80-100 

92 

0-5 

1 

n6 

50-85 

78 

#30 

25-60 

1+9 

#50 

10-30 

16 

#100 

2-10 

1 

#200 

0-3 

1 

Fineness  Modulus  2. 3-3.1  2.62 


TIT.  Specific  Gravity  of  Materials 


Coarse  Aggregate  2.90 

Fine  Aggregate 2.65 

PC 3.15 

IV.  Maximum  w/c  Ratio  6.0* 


*This  is  the  maximum  ratio  allowed  by  specification.  The 
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Table  A 6 (cont'd) 


V.  Class  B Concrete  Weight  Proportions 
1 - 1.65  - 2.96* 

*For  one  pound  of  Portland  Cement,  one  must  have  I.65  pounds  of 
Fine  Aggregate  and  2.96  pounds  of  Coarse  Aggregate. 


Table  AT.  Open-graded  Mix  #1 


I.  Materials 

Crushed  Stone  - Mount  Hope  Materials  (Granite) 
Asphalt  Cement  - AC-20 


II.  Aggregate  Gradation 


Sieve 

Specs . 

Actual 

Size 

(fj  uas sing) 

( % passing) 

2" 

100 

100 

3 A" 

70-95 

80 

3/8" 

40-65 

50 

A 

20-35 

25 

#50 

0-6 

3 

#200 

0-2 

1 

III.  Specific  Gravity 

Aggregate  - 2.68 
Asphalt  Cement  - 1.01 


TV.  Density 

Approximately  131  pcf 

V.  Cylinder  Sizes 

Sample  A 6"  (diam. ) x approximately  8" 
Sample  B 6"  (diam. ) x approximately  8" 

VI . Asphalt  Content 


Unstabilized  Open-graded  Mix  #1. 


% Passing 


Sieve 

Probe 

Guarded 

Size 

Specs . 

Test 

Hot  Plate 

2" 

100 

100 

100 

3/4" 

70-95 

97 

73 

3/8" 

40-65 

68 

43 

#4 

20-35 

35 

20 

#50 

0-6 

5 

3 

#200 

0-2 

3 

0 

Table  A8.  Open-graded  Mix  ti3. 


r 


Table  A9.  Mineralogical  Content  of  Pennington  Traprock 
and  Mt.  Hope  Granite. 


A.  Pennington  Tranrock  from  the  Pound  Brook  Quarry  is  quite  uniform  in 
mineralogy  and  texture.  Mineral  percentage  varies  a little  as  do  the 
effects  of  weathering  (upper  level)  and  hydrothermal  alteration. 

Essentially  the  rock  is  composed  of  the  following: 


1. 

Labradorite  feldspar 

i*o£ 

+ 

2. 

Pyroxene  (mainly  augite) 

1*03 

+ 

3. 

Glass 

5-10 

+ 

U. 

Magnetite 

0-5 

+ 

5. 

Trace  accessories 

Tr 

+ 

There  is  usually  some  alteration  which  forms  the  following: 

1.  Labradorite  - epidote-zoisite , sericite 

2.  Pyroxene  - chlorite 

3.  Glass  - chlorite,  albite 


Fresh  rock  is  gray.  Altered  rock  is  greenish  gray  to  dull  foliage 
green. 

Texture  is  ophitic,  feldspar  crystals  contained  within  pyroxene, 
typical  of  basalt. 

B.  Mt . Hope  Granite 

Alaskite  and/or  Hornblende  Granite  from  Mt.  Hope  is  also  a fairly 
uniform  rock  composed  of : 


1. 

Microperthite  (microcline  with  exsolved  albite) 

1*53 

2. 

Quartz 

25  % 

3. 

Oligoclase 

20 % 

1*. 

Hornblende 

0-10 

5. 

Pyroxene 

6. 

Magnetite 

0-5 

T. 

Misc.  Trace  accessories 

Tr 

The  granitic  rocks  are  quite  uniform,  equigranular , medium  grained 
and  gneissoid. 


APPENDIX  B 


Description  of  Testing 


This  testing  consisted  of  determining  the  thermal  conductivity 
using  the  guarded  hot  plate  method  under  four  different  conditions.  The 
test  conditions  were  as  follows:  Less  than  optimum  moisture  content  at 
95  percent  maximum  density  (a)  above  freezing,  (b)  below  freezing;  and 
field  capacity*  at  95  percent  maximum  density,  (c)  above  freezing,  and 
(d)  below  freezing.  The  physical  properties  (supplied  by  NJDOT)  of  the 
thermal  conductivity  test  samples  are  shown  in  Table  Bl. 


Table  Bl.  Physical  Properties  of  Samples 


Actual  % Passing 


Base 

Course 

Subbase 

Course 

Size 

Ogdensburg 

Gneiss 

Hamilton 

Lakes 

Hamburg 

Shale 

Pennington 

Traprock 

4" 



100 

___ 

. 

o" 

100 

— 

100 

100 

3/4" 

86.9 

98.4 

93.2 

83.1 

#4 

48.0 

74.4 

69.1 

55.1 

9.5 

11.9 

7.2 

l6.8 

#200 

3.8 

1.4 

2.8 

8.9 

Opt i mum 

Moisture 

{%) 

8.1 

7.6 

9.5 

11.1 

Max . Dry 

Density 
(pcf ) 

135.6 

123.6 

127.8 

135.9 

Minimum 

Porosity 

0.210 

0.225 

0.225 

0.244 

*Field  capacity  here  is  the  amount  of  water  a 3-inch  thick  specimen  will 
hold  without  further  drainage. 


Apparatus 

The  apparatus  used  in  testing  the  samples  was  a guarded  hot  plate 
testing  machine,  conforming  to  ASTM  designation  C1TT-71,  ASTM  1973. 
Figures  Bl  and  B2  show  the  physical  apparatus  and  the  test  stack  in- 
stalled, ready  to  be  insulated  and  closed  up  for  test.  Figure  B3  shows 
a cross  section  of  the  stack  assembly  and  the  thermocouple  locations  in 
the  testing  plate  surfaces,  when  assembled  for  the  simultaneous  testing 
of  two  specimens. 
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Figure  B-2.  Test  stack,  installed,  prior  to  insulating 


The  thermal  conductivity  tests  were  performed  in  the  guarded  hot 
plate  apparatus  (Figure  Bl)  in  accordance  with  ASTM  procedures,  except 
for  the  following  modifications.  Instead  of  using  two  identical  samples 
and  evenly  splitting  the  total  heat  flow,  (Q^) , from  the  center  guarded 
plate  through  the  sample,  a standard  gum  rubber  sample  one-inch  thick, 
of  known  thermal  properties,  was  used  to  replace  the  top  soil  sample 
(See  Figure  BM . Then  instead  of  dividing  the  total  heat  flow  (fy.)  by 
two,  the  computed  heat  flow  (Q-,  ) through  the  gum  rubber  was  subtracted 
from  Q-t,  giving  the  heat  flow  (Q^)  through  the  soil  sample.  This  Q2  was 
used  to  compute  the  thermal  conductivity  of  the  soil  sample.  This 
modification  eliminates  the  convective  currents  that  occur  in  the  pore 
water  of  the  top  soil  sample. 


COLD 


□ 

l”  (25.4  mm)  Gum  Rubber  Sample 

HOT 

2"  (50.8  mm)  Test  Sample 

COLD 


Figure  B-U.  Test  Situation. 


Sample  Preparation  and  Testing 

A representative  soil  sample  of  approximately  100  lb  was  mixed  at 
the  desired  water  content  (taken  from  the  compaction  curve)  and  allowed 
to  equalize  for  2h  hours.  Then  a selected  weight  of  wet  soil  was  com- 
pacted to  a pre-determined  volume  in  three  one-inch  thick  layers  in  a 
20"  x 20"  x 3"  deep  mold,  using  an  appropriate  uniform  applied  pressure 
over  the  20"  x 20"  area  of  sample  (see  Figure  B-5).  Prior  to  placing 
the  sample  in  the  Guarded  Hot  Plate  Apparatus  (GKP),  the  sample  was 
frozen  by  first  tempering  it  at  h0°F  for  several  hours,  then  fast  freezing 
it  at  -10°F.  The  frozen  sample  was  placed  in  the  GHP  apparatus  and 
thermal  conductivity  measurements  were  taken  at  approximately  23°F.  The 
sample  was  then  thawed  while  still  in  the  GHP  apparatus  and  thermal 
conductivity  measurements  were  made  at  +1*0°F.  The  sample  was  removed 
from  the  GHP  apparatus  and  brought  to  field  capacity  water  content  by 
supplying  the  bottom  of  the  sample  with  water.  The  wetted  sample  was 
then  frozen  as  before  and  the  same  testing  procedure  followed  as  for  the 
low  water  content  samples. 

The  Computation  of  Thermal  Conductivity,  K 

The  following  is  a short  description  of  what  is  actually  measured 
and  what  is  computed,  using  the  guarded  hot  plate  method  of  measuring 
thermal  conductivity. 
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1 . Known : 

a.  The  test  area  (A)  is  1 ft^  (Figure  B3). 

b.  The  thermal  conductivity  of  the  gum  rubber  sample  (Kgr) 

in  Btu/hr  ft"  °F/in).  This  calibration  versus  temperature 
was  supplied  with  the  gum  rubber  from  the  National  Bureau 
of  Standards . 


2.  Measured: 


a.  The  thickness  in  inches  of  both  the  gum  rubber  (X~R)  and 
the  test  sample  (X  ) , prior  to  starting  the  test. 


b.  The  average  temperature  (°F)  of  each  of  the  four  testing 


plates,  T , T?,  T 
thermocouples  on 


then  averaged  together. 


c.  The  main  heater  voltage  (E)  and  current  (l).  This  is  the 
guarded  heater  within  the  center  of  the  test  stack. 
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3.  Computed : 

a.  The  average  temperature  of  the  gum  rubber  sample: 

TGB  ■ <T1  4 V * 2 

b.  The  average  temperature  of  the  test  sample: 

Tg  = (T3  + Tu)  ♦ 2 

c.  The  differential  temperatures  across  the  Gum  Rubber  and 
the  test  sample: 


d.  The  temperature  gradient  of  the  Gum  Rubber  and  the  test 

sample: 


AT 


GR 


GR 


GR 


e.  The  total  heat  flow  (Qt)  in  Btu/hr  through  both  the  gum 
rubber  and  the  test  sample.  The  mean  heater  power  (watts)  is  computed 
by: 


P = El 


then  using  a constant  conversion  of  3.  *+15  Btu/hr/watt 


Qt  = (P)(3.l»15) 


f. 

from  the  NBS 


The  Thermal  Conductivity  of  the  gum  rubber 
calibration  curve  at  Tntl. 


(kGR)  taken 


g .  The  heat  flow  through  the  gum  rubber  ( ) : 


QGR  = kGRiGRA 

h.  The  heat  flow  through  the  test  sample  (Qg) 

Qs  = Qt  ' QGR 

i.  Finally  the  Thermal  Conductivity  of  the  test  sample  (kg) 

ks  = Qs/isA  in  Btu/hrft2°F/in- 
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Results 


The  results  of  the  thermal  conductivity  tests  are  shown  in  Table  B2 
and  Figure  B6. 

In  all  cases  the  wetted  samples  had  higher  thermal  conductivities. 
Also,  three  of  the  soils  developed  lower  conductivities  when  frozen  than 
when  thawed.  These  results  are  consistent  with  previous  experience, 
which  has  shown  that  thermal  conductivity  rises  generally  with  water 
content.  Samples  having  very  few  fines  will  have  lower  conductivities 
when  frozen  because  there  is  little  unfrozen  water  to  wet  the  mineral 
grains  and  thereby  enhance  the  heat  transfer  from  ice  to  soil  particles. 

In  contrast,  the  Pennington  Traprock  has  considerably  more  fines 
than  the  other  three  soils , and  predictably  a greater  amount  of  unfrozen 
water  at  subfreezing  temperatures.  For  this  soil  the  frozen  condition 
has  a higher  thermal  conductivity  than  does  the  unfrozen  condition. 


W/m  K p«nmngton 


Ttmperatur* 


Figure  B-6.  Thermal  conductivity  of  base  course  and  subbase 

course  materials  from  New  Jersey,  guarded  hot  plate 
method . 
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*These  samples  densified  upon  thawing,  resulting  in  a higher  degree  of  saturation  with  no  change  in  water 
content . 
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Background 


The  thermal  conductivity  of  a soil  is  one  of  the  major  factors  in- 
fluencing the  rate  at  which  freezing  temperatures  penetrate  into  ground 
under  natural  conditions.  Of  particular  interest  is  the  conductivity 
in  the  frozen  zone  adjacent  to  the  location  of  freezing.  Temperatures  in 
this  zone  will  generally  he  in  the  vicinity  of  0°C. 

Kersten  (19^9)  reported  thermal  conductivities  of  various  soils  for 
temperatures  of  25°F  and  UO°F . He  used  a steady-state  method  which  re- 
quired a temperature  differential  of  approximately  10°F.  While  repre- 
senting a considerable  improvement  over  the  Uo°F  differential,  formerly 
required  by  the  guarded  hot-plate  apparatus,  a 10  degree  differential  still 
results  in  errors  caused  by  the  migration  of  moisture  or  vapor  during  the 
many  hours  required  to  reach  a steady  condition.  Latent  heat  effects  at 
temperatures  just  below  freezing  cause  further  difficulties  with  fine- 
grained soils,  in  which  the  amount  of  frozen  water  may  change  appreciably 
with  an  increase  in  temperature. 

Over  the  past  two  decades  a transient  method  of  heating  has  been 
developed  that  reduces  the  intensity  of  these  effects  by  requiring  tempera- 
ture differentials  of  no  more  than  several  degrees..  This  is  the  probe,  or 
modified  line  heat-source,  method.  Basically,  the  method  utilizes  a 
constant  rate  of  heat  production  along  a line  source  to  produce  an  expanding 
cylindrical  temperature  field  in  the  adjoining  material.  The  rate  of 
increase  of  temperature  at  any  point  in  the  material  depends  in  general  on 
its  distance  from  the  source  and  the  diffusivity  of  the  material.  However, 
after  a sufficient  interval  of  time  the  temperature  rise  varies  only 
with  the  conductivity  of  the  material,  and  in  a simple  logarithmic 
manner;  for  points  near  the  heat  source,  corresponding  to  a thin  probe, 
this  period  can  be  reduced  to  15  minutes  or  less  for  most  materials. 
Migration  of  internal  fluids  is  thereby  greatly  minimized. 

Because  of  the  inherent  simplicity  of  the  method,  the  relatively 
small  temperature  differentials,  and  the  short  heating  periods  required, 
the  probe  method  is  well  suited  to  the  measurement  of  the  thermal  conduc- 
tivity of  soils  and  other  materials  near  freezing  temperatures.  Lachen- 
bruch  (1957)  reported  an  investigation  on  frozen  soils  in  place,  in 
which  the  probe  method  was  used.  Also,  a series  of  laboratory  probe 
tests  on  frozen  soils  was  performed  under  contract  for  CRREL  (Wechsler, 
19^5).  Since  that  time,  CRREL  has  been  engaged  in  a continuing  laboratory 
program  with  the  following  objectives:  (l)  to  develop  a simple, 
reliable  apparatus  and  technique  for  performing  such  tests,  (2)  to 
determine  any  special  difficulties  which  may  be  encountered  at  temperatures 
close  to  freezing,  and  (3)  to  compare  measured  conductivities  of 
selected  soils  with  previously  reported  values.  A general  objective  of 
the  program  has  been  to  seek  a physical  explanation  for  the  observed 
thermal  conductivities.  Several  reoorts  are  available  describing  results 
of  this  program  (McGaw  1967,  1968,  1969,  1971,  197^). 
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Theory  of  the  Probe  Method 

The  earliest  mention  of  a line  heat-source  method  of  measuring  con- 
ductivity has  been  attributed  to  Schleiermacher  (l888).  In  this  century 
it  was  suggested  by  Stalhane  and  Pyk  (1931 ) as  a means  of  measuring  the 
thermal  conductivity  of  a liquid.  The  first  practical  use  of  the  method 
is  credited  to  Van  der  Held  and  Van  Drunen  (19^9)  who  used  a heated  wire, 
also  in  connection  with  liquids.  Both  the  heated-wire  method  and  the 
probe  method  are  based  on  the  simple  line  heat-source  theory. 


Line-Source  Theory.  Beginning  from  an  initial  isothermal  condition, 
and  with  a constant  heating  rate  (q)  per  unit  length  of  source,  the  simple 
line  heat-source  theory  results  in  the  following  well-known  expression 
for  the  temperature  rise  at  a distance  r from  the  line  source  during  the 
heating  time  t_: 

2 


Tr  = nS  Ei  {VTT]  ’ 


(1) 


where  a is  the  thermal  diffusivity  of  the  external  medium,  and  K is  its 
thermal  conductivity  (Carslaw  and  Jaeger,  1959)-  The  exponential 
integral  (Ei)  is  tabulated  by  Jahnke  et  al  (i960),  but  is  also  commonly 
approximated  by  the  series  expansion, 


-Ei 


2 


in  Y _ (^)  + i/l*  (f^)  + 


(2) 


where  y = 0.5772.  When  t is  sufficiently  large  (or  r is  sufficiently 
small)  that  the  value  of 


at 


•»1/1*  , 


then 


1 ’ 


(3) 


and  all  but  the  first  two  terms  of  the  expansion  are  considered  to  be 
negligible.  Hence,  equation  (l)  becomes 


T = JL_  (m  J*«L  _ Y)  . 

r UttK  V 2 
r 


or 


t = In  t + (in  — - y), 
r l*irK  2 

r 


(1+a) 


(Ub) 
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With  a and  £ unchanging,  the  expression  in  parentheses  is  a constant. 
Consequently,  in  the  simple  theory  the  temperature  rise  plotted  against 
the  logarithm  of  time  will  be  a straight  line  when  the  criterion  given 
in  equation  (3)  holds. 


In  practice,  the  temperature  rise  (AT)  between  any  two  large  succes- 
sive times  (t  and  t^)  is  utilized,  in  which  case  equation  (Ub)  becomes 


(5) 


the  value  of  the  conductivity  being  given  by  the  basic  line-source 
equation: 


K = 


(6) 


Probe  Theory.  Hooper  and  Lepper  (1950)  were  the  first  to  use  a probe , 
or  portable  heating  element,  for  measuring  the  conductivity  of  soils  in 
their  natural  condition.  They  were  followed  shortly  after  by  de  Vries 
(1952),  Lentz  (1952),  Hooper  and  Chang  (1953),  and  others.  These  in- 
vestigators found  that  definite  precautions  must  be  taken  in  the  design 
and  operation  of  a probe  if  it  is  to  approximate  a simple  line  source. 

The  dimensions  of  the  probe,  its  thermal  properties,  and  the  contact 
resistance  between  the  probe  and  the  test  material  all  must  be  taken  into 
account,  either  directly  or  indirectly.  Consequently,  beginning  with 
D'Eustachio  and  Schreiner  (1952),  an  extensive  series  of  analyses  of  the 
possible  errors  in  the  probe  method  was  carried  on  by  Van  der  Held  et  al 
(1953),  Blackwell  (1951*),  Vos  (1955),  3uettner  (1955),  Jaeger  (1956), 
Blackwell  (1956),  Jaeger  (1958),  de  Vries  .and  Feck  (1958),  and  others. 

A summary  of  these  studies  is  given  by  Wechsler  (1965). 

Typical  solutions  for  the  temperature  rise  at  the  inner  radius  of  a 
hollow  metal  probe,  heated  along  its  outer  radius  b,  have  been  given  by 
31ackvell  (195*0.  If  the  probe  ^as  a heat  capacity  M.C^  per  unit  of 
length  and  a contact  resistance  — per  unit  of  area  between  its  outer 
surface  and  the  external  medium,  the  large-time  solution  is: 


[ In  Ut  - y 


2K 

bK 


+ 


_1 

2t 


(7) 


[in  Ut  - y + 1 - j (In  Ut  - y + + 0 (-|-)]  , 

2 

where  t = — and  8 = (yy,"—  , pc  being  the  volumetric  heat  capacity 
b2  M1C1 

of  the  external  medium. 


r 


Jaeger  (1956)  found  an  identical  solution  for  a solid  cylindrical 
probe  under  the  same  conditions.  Consequently,  whether  a probe  with  an 
internal  heater  is  considered  to  be  hollow  or  solid  makes  no  change  in 
the  large-time  solution.  Moreover,  a related  solution  by  Blackwell  (195*0 , 
in  which  he  takes  into  account  the  finite  conductivity  of  the  probe 
material,  is  also  identical  to  equation  (7)  except  for  the  addition  of 
one  other  term: 


57  t-rf<4i*42» 
0 


(8) 


here,  and  A0  are  functions  of  the  inner  and  outer  radii,  a and  b. 
According  to  Blackwell,  this  term  is  usually  insignificant  in  magnitude. 

It  therefore  appears  that  the  finite  conductivity  of  a probe  also  pro- 
duces no  significant  change  in  the  large-time  solution,  so  long  as  the 
probe  material  is  a good  conductor  in  relation  to  the  external  medium. 

An  axial  flow  of  heat,  brought  about  by  temperature  gradients  along 
the  probe  or  by  conduction  along  the  electrical  or  thermocouple  leads,  may 
introduce  serious  distortions  from  equation  17)  by  reducing  the  effective 
heat  input  by  an  unknown  (and  perhaps  varying)  amount.  Lentz  (1952) 
and  Hooper  and  Chang  (1953)  took  note  of  conduction  along  the  leads, 
while  Blackwell  (1956)  analyzed  the  axial  flow  error  in  the  probe  itself. 
Blackwell  found  that  temperature  gradients  along  a carefully  constructed 
probe  of  good  conducting  material  should  be  negligible  if  the  length  of 
the  probe  is  at  least  25  times  the  diameter.  In  practice,  length-to- 
diameter  ratios  of  100  or  more  are  usually  employed,  so  that  axial  heat 
losses  are  generally  assumed  to  be  negligible.  On  the  other  hand,  loss 
of  heat  from  the  system  through  the  leads  is  a distinct  possibility 
whenever  the  ends  of  the  probe  are  at  a higher  temperature  than  the 
surroundings.  This  latter  factor  appears  to  have  been  overlooked  in  the 
more  recent  literature. 

Assuming  that  adequate  measures  have  been  taken  to  control  axial 
heat  losses,  equation  (7)  may  be  accepted  as  the  basic  probe  equation, 
in  which  contact  resistance  and  the  thermal  properties  of  the  probe 
have  been  taken  into  account.  As  in  the  line  heat-source  analysis, 
incorporating  the  criterion  of  equation  (3)  will  eliminate  terms  pro- 
portional to  powers  of  (y)  • That  is,  for  times  sufficiently  large 
that  x >>  1/**,  the  temperature  rise  at  the  center  of  the  probe  is  given 
t>y 


t » -i-  rln  hsL  . Y + aL] 
LttK  1 v2  Y bH J 
b 


or 


Uwk  [ln  + (ln  rl- Y + bl1 

D 


(9a) 


(9b) 


65 


i 


r 


Since  b and  y_  are  constants,  the  expression  in  parentheses  in  equation 
(9b)  is  a constant  when  a,  PC,  and  (^)  are  unchanging  during  a test 

period.  In  this  case  the  temperature  rise  between  successive  large 
times  reduces  to  equation  (5)»  the  value  of  the  conductivity  being 
given  by  equation  (6)  as  before.  Thus,  in  a well-designed  test  the 
basic  line-source  equations  also  hold  for  the  probe  method. 

Equation  (9b)  may  also  be  written  in  the  following  form,  which  is 
sometimes  convenient  when  measuring  the  conductivity  of  a frozen  soil 
near  the  freezing  point: 

T = [ln  t + (ln  72  " y)  + ln  £ + bH]  ‘ (9c) 

b 

Here,  the  term  in  parentheses  is  an  absolute  constant,  while  the  last 
two  terms  may  each  be  variable.  There  is  some  danger  that  a progressive 
change  in  either  term  during  a test  will  give  an  erroneous  slope  to  the 
graph  of  T_ versus  ln  t while  not  noticeably  affecting  the  linearity. 
Consequently,  special  care  must  always  be  used  in  the  interpretation  of 
test  data  for  frozen  soils . 

Apparatus 

The  probe  used  at  CRREL  was  constructed  by  Arthur  D.  Little,  Inc. 
of  Cambridge,  Massachusetts.  It  is  similar  in  construction  to  that 
used  by  D'Eustachio  and  Schreiner  (1952),  although  somewhat  smaller  in 
diameter  and  approximately  twice  as  long.  A stainless  steel  sheath 
of  0.065  in.  outside  diameter  contains  a heating  coil  of  constantan 
wire  (66.0  ohms/in.)  within  which  is  placed  a chromel/constantan  thermo- 
couple. The  thermocouple  lies  on  the  axis  of  the  probe  and  at  its  mid- 
length. A molded  plastic  block  at  one  end  supports  the  thermocouple 
and  power  leads.  The  free  length  of  the  probe  below  the  block  is  8.00  in., 
resulting  in  a length/diameter  ration  of  125.  A schematic  of  the  probe 
construction  is  shown  in  Figure  Cl. 

The  probe  circuitry  consists  of  a heater  circuit  and  a thermocouple 
circuit.  The  heater  circuit  (Figure  C2a)  contains  the  probe  coil  of 
ohms,  a switch,  a DC  milliameter  sensitive  to  0.001  ma,  and  a 
regulated  DC  voltage  power  supply.  The  heat  input  to  the  coil  is  read 
as  12R. 

The  probe  thermocouple  circuit  (Figure  C2b)  contains  the  chromel/ 
constantan  probe  junction,  a reference  junction  placed  in  an  ice  bath, 
a bias  voltage  (accuracy  0.1?!),  and  indicating  amplifier,  and  a milli- 
volt recorder.  At  any  time  the  probe  junction  temperature  can  be  read 
directly  as  the  algebraic  sum  of  the  bias  voltage  (+)  and  the  scale 
voltage  ( + ) as  shown  either  on  the  amplifier  indicator  or  the  recorder. 
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The  recorder  is  operated  as  zero-left  with  +20  microvolts  full  scale 
(0.5F).  Prior  to  a test  the  bias  voltage  is  adjusted  to  maintain  the 
trace  on  the  recorder  scale  from  approximately  10  seconds  of  heating.  In 
most  cases  linearity  of  the  trace  will  have  been  established  before  60 
seconds  of  heating. 

Standard  test  specimens  are  three  in.  (76  cm)  in  diameter  and 
approximately  8-1/2  in.  (63  cm)  long.  They  are  molded  in  an  aluminum 
cylinder  12  in.  (30  cm)  long  having  a wall  thickness  of  1/4  in.  (0.63  cm) 
and  a rigid  insulating  bottom  (Figure  C3).  This  container  fits  snugly 
into  an  aluminum  cooling  well  bathed  by  circulating  pre-cooled  ethylene 
glycol  (Figure  C4).  A sealed  air  gap  surrounding  the  circulating  chamber 
provides  insulation  from  the  room  temperature,  which  is  maintained  at 
+15°C  (59°P).  One- inch  thick  foam  insulation  is  placed  beneath  and 
above  the  cooling  well  assembly.  Insulation  is  also  applied  to  the  top 
of  the  specimen  container  after  it  has  been  placed  in  the  well  and  the 
probe  inserted. 

Temperature  regulation  of  the  ethylene  glycol  refrigerant  is  ac- 
complished in  an  external  cooling  bath,  in  normal  operation.  In  the 
case  of  the  New  Jersey  specimens,  which  were  6 in.  in  diameter  because 
of  the  large  aggregate,  a somewhat  modified  procedure  was  devised.  In 
place  of  the  standard  container,  a waterproof  plastic  cylinder  was  fitted 
snugly  over  the  specimen,  enclosing  the  sides  and  the  bottom  of  the 
specimen.  This  assembly  was  then  suspended  directly  in  the  ethylene 
glycol  refrigerant,  thereby  cooling  the  specimen  to  the  regulated  bath 
temperature. 

After  insertion  of  the  probe,  a 1-in.  thick  disc  of  styrofoam  in- 
sulation was  fitted  over  the  specimen  to  protect  it  from  the  room  air, 
which  is  maintained  at  40°F.  Probe  heating  was  then  accomplished  in 
the  normal  manner,  and  thermal  conductivity  was  calculated  from  the 
data  as  described  in  the  next  section. 

Evaluation/Calibration  Procedure 

Evaluation  of  the  probe  apparatus  is  accomplished  at  a temperature 
of  +5.0°C  using  three  materials  whose  conductivities  are  well  known. 

These  materials  are  Ottawa  Standard  Sand  (20-30  mesh),  both  dry  and 
saturated,  and  distilled  water. 

Typical  results  obtained  for  water,  dry  sand,  and  saturated  sand 
at  +5°C  are  given  in  Figures  C5,  C6,  and  C7.  From  the  value  of  the 
measured  slope  (AT  vs_  In  t),  thermal  conductivity  is  calculated  from  the 
following  simple  relation: 


K = 


(10) 


68 


cylinoer 

Aluminum 


lOOOcm3 

SOIL 

SPECIMEN 


aluminum 

cylinoer 


MONITORING 

thermocouple 


POROSITY:  0.355 
WATER  CONTENT:  21.4% 


0.1  0 2 0 4 0 6 0 8 I 2 4 6 8 10 

TIME , minutt* 

Figure  C-T.  Initial  Trial:  Saturated  Ottawa  Standard  Sand. 

where  I is  the  measured  current  in  amperes , AT  is  the  temperature 
difference  along  the  linear  portion  of  the  graph  (in  microvolts)  taken 
over  a single  decade  of  time  (0.5  min.  to  5.0  min.;  1.0  min.  to  10.0  min., 
etc.),  and  C is  a constant.  Equation  (10)  is  derived  from  the  probe 
constants  and  equation  (6)  as  follows: 

A-  K = TT^at  ln  ; (6) 

2 

3.  q = I H,  per  unit  length  of  probe,  where  R = 26.0  ohns/cm. 

= (26.0  t2  x (0.2389  p-~/^c  ) 

cm  watt 

2 

= 6.21  I cal/cm-sec; 
t t t 

C.  In  = 2.303  log. - = 2.303  (for  -r-  = 10) 

1 XU  rl  Z1 

D.  Combining: 

v — 6.21  /p  _ 1 • 13 1 j2  cal 

^'l^  1 (2.303)  - I uv.cm_sec 

E.  58.5  pv/°C  ( chromel/constantan ) § 0°C 

F K . (1.137  (58.5)  2 . 66^5  ,2  cal  , . 

K ATX  1 ATj_  1 °C-cm-sec  * U0  > 


The  tests  for  water  typically  reach  good  linearity  prior  to  30  seconds 
of  heating.  With  the  lowest  level  of  heat  input  (1.5  mw/in. ) convective 
disturbances  do  not  set  in  until  after  5 minutes  of  heating;  with  higher 
heat  inputs,  convection  begins  after  about  2 minutes.  Both  heating  and 
cooling  curves  give  similar  results.  In  three  tests  the  following  con- 
ductivity values  were  obtained:  1.33,  1.37,  1.3^  (avg.  1.3^*7,  heating); 
1.25,  1-37,  1.3^  (avg.  1.320,  cooling).  The  value  at  +5°C  is  reported 
to  be  1.3^  mcal/°C-cm-sec . (Dorsey,  19*+0,  Table  130  (BN)). 

There  is  no  problem  from  convection  in  the  dry  sand.  Linearity 
is  established  prior  to  12  seconds  on  heating  and  30  seconds  on  cooling. 
In  five  trials  the  following  conductivity  values  were  obtained:  0.71, 
0.72,  0.77,  0.73,  0.72  (avg.  0.730,  heating);  0.75,  0.69,  0.67,  0.72, 
0.72  (avg.  0.710,  cooling).  Kersten  (l9*+9)  found  0.77  mcal/°C-cm-sec 
for  the  same  density  at  U0°F  by  a steady-state  method;  de  Vries  and  Peck 
(1958)  found  0.73  at  +17. °C  by  a transient  method  using  a heated  wire  in 
a glass  capillary. 

The  tests  using  water  and  dry  sand  show  that  adequate  precision  in 
the  computation  of  conductivity  is  obtained  when  the  slope  of  the 
linear  portion  (AT^),  over  one  decade  of  time  is  at  least  1 uv.  In  the 
course  of  the  saturated  sand  trials  it  was  found  that  from  four  to  five 
times  the  heat  input  was  required  to  obtain  a value  of  AT^  equal  to 
1 uv.  Woodside  and  Messmer  (1961)  found  the  conductivity  of  saturated 
20-30  Ottawa  sand  to  be  about  8.0  mcal/°C-cm-sec , or  10  times  that  of 
dry  sand.  Trials  with  saturated  Ottawa  sand  give  a value  of  8.22  meal/ 
°C-cm-sec . 
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MODIFIED  llFRGGREY  SOLUTION  FOR  MULTILAYERED  SYSTEMS  (GWA  VERSION) 
G.  AITKEY/R.  HE  R3,  MARCH, 1966 


* LATEST  UPDATE  4/9/63,  & WA 

* ’H  I S VE  HS  I •••.  Pt;  'V  I DE  S F 'R 

* 1.  AL3HA  HEADS*  AND  SJMMARY  DATA 

* 2.  TYSE  IN  OF  3ASIC  DATA 

* 3.  KICK  OUT  FDR  DEPTH  CHANGE  LESS  THAN  .03  FT 

* OPERATING  INSTRUCTIONS 

* A.  SET  F L A 3 1 = 1 . FOR  FROST  DEPTH  SOLUTION,  F.A51=3. 

* F ~R  'HAW  DC  P M S LU  TI  "•  . 

* «.  SET  FLAG2=1.  for  FILE  INPUT  ; FLAG2=0.  for  terminal 

* INPUT. 

REAL  KIND 

C HA  SAC  E R»  3 A ME 

101  PRINT,"  INPJT  F L A 3 1 = 1 . IF  LOOKING  FOR  FROST  DEPTH’* 

I NPU  T , F L A3  1 

PRINT, "ARE  YOU  READING  FROM  A FILE-  TYPE  '1.*  IF  YES'* 

READC0,1)  FLAG2 
IF  (FLAG2.E3.1  .)  30  TO  111 
GO  TO  112 

111  PRINT, "IN’UT  AN  8-LETTER  FILENAME" 
u E A D C 0 , 1 ) "A  M E 

OPEN  FILE  1 , N A M E 

112  A1  =.i)7052337R4 
A2=. 04228201 23 
A 3=. 0092735272 
A 4 = . 0001  5201  43 
A 5-  .0002755672 
A6= .0000450638 

109  IF  ( FLAG2. EQ . 1 .)  30  TO  104 
PRI A • 108 

108  FORMATCIh  ,58  m TYPE  OtGREF-DAYS  < DE G R E E ) , N - F A C T 0 R ( E N ) , L E N G T H 0 F FREEZING'  ' 
K52H  OR  THAWING  SEASON  (T),MEAN  ANNUAL  TEMPERATURE  (ATM)) 

READ(0,1)  DEGREE, EN, T, ATM 
G T " 1 1 7 

1 J 4 READ(1,1)  DEGREC,EN,T,.*TM 
1 F 0 R iM  A T ( V ) 

117  CONTINUE 
DE  P'H  ~=0  . 

DEPTH® J. 

SUMD=0. 

SUMFD=0. 

SUMC  D = 0. 

SUM  I ND  = 0. 

RES=0. 

S UMR  ES  *0  . 

«*• =10 
N *0 
S I = 0 . 

105  DFGREE»AgS( DEGREE) 

V A 6S  ( A ’M  - 32  . ) 

VS*EN* DEGREE /T 
DEGREE*0E3REF.  » E N 
ALP=VO/VS 
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91 D 2 

913 

914 
91  5 
916 
91  7 
91  8 
91  9 

920 

921 

927 

928  1 1 3 

929 

930  114 

931  115 

932 
940  4 
94  5 * 
950  6 
960 
970 
930 
935  * 
993  7 
1000 
1010 
1020  20 
1030 
1040 

1 050  5 


1060 

21  7 

PRINT  208 

1070 

2 08 

FTRMA  T(  •!  NPU  1 C,  THERM,  FUSI”*.*) 

1071 

* 

C - HEAT  CA=>ACITY 

1072 

* 

THERM  - THERMAL  CONDUCTIVITY 

1073 

* 

FUSION  - HEAT  OF  FUSION 

1080 

209 

HEAD(0,9)  C,  ’ H E R M , F U S K' 

1090 

9 

FORMAT  (V) 

1100 

GO  TO  21 

1110 

220 

RE  A D ( 1 , 9)  C , T HE RM, FUS  ION 

1120 

21 

A * 'HE  R-M  /C 

1130 

SUM  D * SU  MO  HD 

1140 

depth=sumd 

1150 

10 

SUMFD=SUMFD‘(FUSI3N*D) 

1160 

BAPL  = SUMF  D/SUNP 

1170 

SUMCD*SUMCDH (C‘D) 

1180 

oarc=sum: d/sumd 

1190 

IF(BARL)S7,37,36 

1200 

36 

U*VS‘(mARC/BARL) 

1210 

G A R B *1 ,77245/U 

1220 

R * . 1 5 

1230 

30 

ERF*1.-1./((1.‘A1*R‘A2‘R“2*AS‘R“3‘A4»R“4‘A5»R“5 

1240 

K »A6‘R“6> “ 16) 

1250 

EXPF  = EXP(-(R“2)) 

1260 

ERFC  = 1 .-E  RF 

1270 

CHECK*EX»F/ERF-((4LP)‘EX°F/ERFC) 

1280 

CHECKAsR-*  GA  R8 

1290 

IF(FLAG)  1002,1  002,1  00.3 

1300 

1 U02 

FL AG«1 00. 

FL  AG  = -1 00 . 

FL AGAa-1 0J  . ‘ 

IF  (FLAG.?. EU. 1.)  30  TO  113 

PHI  5 I TP  J • KIf.O  (4  *EAL  NUMBER  - LESS  H4  ' 2.  MEANS  'HE  SOIL" 

PRINT,"  IS  COARSE-GRAINED,  GREATER  THAN  OR  EQUAL  TO  ?.  BUT" 

PRINT,”  LESS  THAN  10.  MFANS  THE  SOIL  IS  FINE-GRAINED,  AND" 

PRINT,"  GREATER  THAN  OR  EQUAL  TO  10.  MEANS  YOU  V I L L HAVE  TO" 

PRl?.,"  SPECIFY  ADOI'I—AL  THERMAL  parame -er-s>  ,GAMMA  (DENSITY)," 

PRINT,”  WC  (HATER  CONTENT,  PERCENT  ON  A GRAVIMETRIC  BASIS),” 

PRINT,”  AND  D (LAYER  THIC<NESS,  IN  FEET)” 

GO  TO  114 

READ  (1,1)  KI  ’ D, GAMMA, WC,D 
GO  TO  115 

READ(0,1)  <IND, GAMMA, WC,D 
CONTINUE 

I F(< IAO-1 0 . ) 4 , 5, 5 
IF  (KIND-2.)  6,7,7 
COURSE-GRA  IN ED 

THERMT  = (.  7*ALOG10(VC)  + .4)  * ( 1 0 . * * ( .01*  GAMMA) ) 

■'HERMF  = .076*(1U.**(.013*GAIMMA))4.0  32*(1C.»*(.0146*GAMMA))*UC 
THERMa(THE  RMT  + THFRMF ) /24 . 

GO  TO  20 
F INE-GRA INED 

*HEKMt=(.  9*AL  .’G10(  JC  )-.2)  * (10.*  *(.01*  GAMMA  ) ) 

THERM F = .01*(  10.**(.022*GAMMA)  ) + .085*(  1 0 . * * ( .008*GAMMA))*WC 
THERM=(THERMF4THERMT)/24. 

C«GAMIMA  * ( . 1 7*  (.  75*riC/1  00.  ) ) 

FUSJ' ',=  144. ‘GAMMA ‘(WC/103.) 

GO  TO  21 

IF  (FLAG2.EQ.1.)  GO  TO  2?0 
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1510 
1 320 
1 5 5 J 
1 340 
1333 
1360 
15  70 
1333 
1390 
1 430 
1413 
14  2 3 
1433 
1 440 
1 453 
14  63 
1465 

1473 

1474 

1475 
1 43  J 

1490 

1491 

1492 

1493 

1494 

1495 

1496 
149  7 

1498 

1499 

1500 
1 5 31 
1502 
1533 
1 534 

1535 

1536 

1537 
151  J 
1520 
1530 
1543 
1553 
1 563 
1570 
1530 
1590 
1600 
1610 
1623 
1 630 
1643 
1 650 
1 663 
1673 
1 680 
1690 


IF  (CHE C< -CHECK  A)  1 3 04,1030, 1 nOU 
1004  FL AGAslOO. 

GO  TO  10J1 

1J03  I F ( FLAGA)  1003/  1000/  1 001 

1000  IF(CHFC«-CHCCK4)35,35,34 

1001  IF (CHECK  — CHECKA)  100  7»35#35 

1007  R=R-.0J1 

GO  TO  30 

34  -7  * i<  ♦ . 0 3 1 
60  TO  3C 

35  F=((2.*(R*«2))/J)*».5 
60  TO  38 

37  F=0. 

0*0  . 

33  WES=A8$(D)/THF7M 

T0TAt_R  = SJMRES  + RES/2. 

IF  (PAFl)  40/40,3> 

39  AYERI=(FJSI0N*(A3S(D))/(24.»(F»*2)))*T3TAIR 

60  TO  41 
43  AYERI=0. 

41  SI=SI+AYERI 

N«N  »1 
PR  I NT 

IF  (N.EQ.1)  PRINT  200 

2 JO  F -RF'A  1 ( 1 H /25X/53HMUL  TILAYER  S LUTri  F THE  MODIFIES 

& PE  R6G  R b N E3UAT  I 3 N ) 

IF  (N.E9.1)  PRINT  136 

106  FORMAT (1 H /I 1HALPHA  INPUT) 

IF  ( '..  E 0.  1 ) PRIT-  51  J 

510  FORMAT  (1H  ,//) 

IF  (N.EQ.1)  PRINI  931 

201  FORMAT  (1  H ,11X,3HBRY,4X,SHiJATER,4X,5HL4YE»,4X,4HYeAr,4X,7HTHERNAL, 

S2X/6HLA  E ' *,  7h  FJSI  * ,7H  L A M H DA  ,2  X , 5H  l A Y E «•,  3X  , 4 H C J M.  , 6 X , 5 H L A Y E R,  4 X 
i,4HCUM.,/,10X,7HDENSITY,3H  CONTENT, 3x,4HSIIE,3X,3HCAP4CITY,iX,5HC0 
4NB.,4X,4HHEAT,3X,4H3AR.,11X,4HRES.,3X,4-(RES.,5X,5HIN0EX,3X,5HIN0EX 
i,/,mx,7H(6A'|MA),2X,4H(UC),6X,3H(0),5X,3H(C),5X,7H(THERM),9H  (FJSI 
K ),2x,3H(U)  ,3X,3H(F  ),4X,SH(RES),9h  (SUMN€S),2x,7h(AYERI),3x,4H(SI 

%),/) 

PRINT  11,  N ,3AM M A, 4C/9/C/THERM, FUSION, U,F,RES/SJMRES,4YERI/SI 

11  FORMA T(  3X, I 3, 2X,?F 8.2,2 X,F  7.3,2 X,F6 .2,2 X, F3. 4,2 x,F 6 .0,2 X,F5. 2, 

It  2X,F4.2,2x,Ffc.?,2x,F6.?,2x,F7.0,2x,F7.3) 

I f (NNN) 15 ,1 2, 12 

12  IF(DFGREE-SI)15, 133,14 

14  IF((DE3REE-SI>-n.)100, 100,501 

501  SlIMPES  *S  JM  re  S ♦ R"F  S 

GO  TO  2 

13  C ON  *S  U M D - D 
BEPTH9*DEPTH-B 
S I Al«S I -A  YE  R J 
NNN  s- 1 0 

15  I F(OEGRFE-S I ) 1 7, 1 03,1 6 

16  I F( DEGREE- ( SI »1 3. ) ) ICO, 100, 18 

18  SUMFD*SUMFO-( FUS  I . ' *0) 

SUMCD=SUMCD-(C«B> 

B=0*((0E3REE-SIN1 )/AYERI) 

IF((AHS(SUM0-(CON*0)  ))-.()  3)100, 100, 132 
102  SUMD  = C -.'♦0 

SI *S I M 1 


79 


in: 

171  0 
17*0 
1 730 
1 74  0 
1 750 
1 760 
1773 
1750 
1 790 
1 3 JO 
1510 
13*0 
1330 
1840 
1 350 
1 360 
1370 
1880 
1390 
1900 
1 91  3 
1915 
1 9*0 
1 933 
193* 
19J5 
1939 
1 943 


GO  TO  10 

17  IF<(DE3REE+10.)-SI)1?,100,10C 

100  DEPTHsDEPThO+D 

PRINT  7Q1 .DEGREE 

701  FORMATdH  ,////, 5X.8HINDEX  = , F 7 . 0 . 1 X , 1 3 HD E GR E E DATS  F) 

PRI'T  702. F- 

70*  FORMAT (1  H .4X.11HN  FACTOR  = . F 5 . 2 ) 

IF  CFLAG1  .NE  . 1 . ) 30  TO  33 
PRINT  703. DEPTH 
G"  747 

33  PRINT  203.DCPTH 

*03  FORMAT(H  ,4x,465DE3TH  OF  THAW  PENETRATION  FOR  THIS  SOLUTION 
8 = .F4.1.1X.5HFEET.) 

703  F^RMAMIH  .4x,47hDEP'H  *F  FR’ST  HENE  HATI—  F R THIS  SOLUTION 
8 * .F 4.1. IX, 5HFEE  T . ) 

747  PRINT  704, T 

704  FORMAT (IN  , 4 X , 1 9HL E N G TH  OF  SEASON  * , F 5 . 0, 1 X.5HD4YS.) 

P K I T ” 705 , A ’M 

705  FORMAT  (H  ,4X  ,19-dEAN  ANNJAL  TEMP  * .F5.1.1HF/) 

500  PRINT  *04 

*04  FORMATE  'SOLUT  ION  COMPLETE.  TYPE  "100003."  If  YOJ  DO  NOT  WANT*,/ 
8*  **  C * T I *.U  E A D "99999."  IF  YOU  DO.') 

1999  IF  (FLAG*. EQ . 1 . ) READ(I.I)  TAG 
IF  (FLA32.NE.1 .)  READ(O.I)  TAG 
IF  (TAG. EO. 100000.)  G*  T"  *001 
IF  (T AG .E 3 .99999. ) GO  TO  109 
GO  TO  500 
*001  CONTINUE 
END 


